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ABSTRACT 

A  geologic  map  of  part  of  the  Brazeau  Range  Anticline  in  the  Rocky  Mountain 
Foothills  near  Nordegg,  Alberta  has  been  made  on  the  scale  of  4  inches  to  1  mile. 

The  results  of  the  mapping  are  presented  as  structure  contour  maps  of  the  tops  of 
the  Alexo,  Palliser  and  Banff  Formations. 

In  the  preparation  of  the  structure  contour  maps  a  new  method  for  determining 
the  generatrix  of  a  cylindrical  fold  based  on  a  generalization  of  the  normal  distribution 
onto  a  sphere  has  been  used.  This  method  allows  a  surface  to  be  subdivided  into  domains 
which  have  been  cylindrically  folded.  Orientations  of  the  bedding  surface  are  then 
extrapolated  down  the  generatrix  of  the  fold  from  exposed  areas  of  the  surface  to  give 
information  on  covered  areas  of  the  surface  within  the  same  domain. 

A  kinematic  analysis  of  the  northern  half  of  the  Brazeau  Range  suggested  the 
Range  is  a  northwesterly  plunging  anticline  formed  by  the  northeastward  movement  of 
strata  up  the  Brazeau  Range  Fault  -  a  splay  of  the  Brazeau  Thrust.  The  structures 
observed  in  the  thesis  area  conform  with  this  analysis.  In  the  dynamic  analysis,  the 
form  of  the  thrust  faults  which  develop  in  the  plane  strain  of  a  perfectly  plastic  body 
moving  across  a  perfectly  rough  surface  has  been  calculated.  These  resemble  the 
thrust  faults  which  developed  in  the  Foothills  Belt  of  Central  Alberta. 
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INTRODUCTION 

In  this  thesis  an  attempt  is  made  to  trace  the  steps  involved  in  the  structural 
analysis  of  a  small,  well-exposed  area  of  sedimentary  rocks  where  the  structure  is 
probably  the  result  of  only  one  episode  of  deformation.  The  analysis  is  in  three 
parts,  geometric,  kinematic  and  dynamic. 

The  geometric  analysis  determined  the  orientation  of  the  bedding  planes  to 
the  exclusion  of  other  surfaces  and  lineations  (such  as  slickensides  and  joints)  of 
possible  structural  significance.  The  results  are  given  as  contour  maps  on  parts  of 
three  stratigraphic  surfaces.  These  surfaces  have  been  subdivided  into  domains  - 
portions  of  the  surface  statistically  determined  to  be  cylindrical ly  folded. 

The  kinematic  and  dynamic  analyses  are  based  on  the  assumption  that  the 
geometry  of  the  surfaces,  the  movements  that  placed  the  surfaces  in  their  present 
position,  and  the  stress  field  that  caused  the  movements  all  have  the  same  symmetry. 
This  assumption  leads  to  difficulties,  for  the  geometric  analysis  reveals  a  large 
number  of  small  domains  in  the  structure. 

A  kinematic  analysis  has  been  put  forward  but  rather  than  attempt  the  trivial 
solution  of  the  dynamic  analysis,  the  thesis  begins  work  on  a  more  general  approach 
aimed  at  discovering  not  only  the  orientation  of  the  principal  stresses  but  also  their 
size  and  history.  There  is  a  brief  discussion  of  the  mathematical  models  that  have 
been  proposed  for  deformation  of  the  type  shown  in  the  Brazen* ■  Range,  and  the 
analysis  of  a  new  model  based  on  the  assumption  that  the  rocks  involved  in  the  defor 
motion  behaved  as  a  perfectly  plastic  body  moving  across  a  perfectly  rough,  horiz¬ 
ontal  surface  towards  a  stable  area  has  been  started.  The  analysis  shows  that  the 
type  of  fault  which  would  develop  under  these  circumstances  closely  resembles 

faulting  observed  in  the  Foothills  region. 
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The  eventual  outcome  of  the  analyses  of  mathematical  models  will  be 
predictions  about  the  form  and  location  of  bedding  surfaces  and  faults.  If  these  are 
to  be  compared  with  field  evidence  the  limits  of  accuracy  of  mapping  and  extra-= 
polation  into  covered  areas  should  be  known.  The  thesis  begins  the  investigation 
of  these  limits. 


Summary  of  Previous  Work 

Early  geological  work  in  the  Nordegg  area  is  described  in  detail  by  Sanderson 
(1958)  and  Pitcher  (1958).  The  first  work  carried  out  by  the  Geological  Survey  of 
Canada  in  the  area  involved  the  mapping  of  the  Mesozoic  coals  of  the  Nordegg  and 
Bighorn  basins.  Systematic  mapping  of  the  thesis  area  by  the  G.S.C.  began  in 
1940.  The  map  areas  along  the  Brazeau  Range  were  first  mapped  on  the  scale  of 
2  inches  to  1  mile,  and  then  revised  on  the  scale  of  1  inch  to  1  mile.  The  revised 
maps  of  the  thesis  area  are  Nordegg  (Douglas,  1956),  Chungo  Creek  (Douglas,  1958) 
and  Alexo  and  Saunders  (Erdman,  1950).  The  thesis  area  also  includes  the  extreme 
southwest  corner  of  the  Harlech  map  area  which  the  G.S.C.  has  yet  to  map. 

In  1958,  the  Alberta  Society  of  Petroleum  Geologists  held  their  eighth  annual 
field  congress  in  the  Nordegg  area.  The  guidebook  contains  a  4  inches  to  1  mile  map 
compiled  from  published  and  unpublished  mapping  of  about  1,000  square  miles  centred 
on  Nordegg  townsite.  With  a  section  and  notes  by  A.N.  Thomas,  it  gives  a  general 
picture  of  the  area.  Shaw  (1963)  contained  a  section  through  the  area.  An  index 
map  of  the  area  showing  G.S.C.  mapping  forms  Fig.  1  . 
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Location  and  Accessibility 

The  area  mapped  for  this  thesis  includes  the  southeast  end  of  Shunda  Mountain, 
the  whole  of  Coliseum  Mountain  and  continues  southeastwards  along  the  Brazeau 
Range  to  take  in  part  of  the  unnamed  hills  which  overlook  Shunda  Gap  from  the 
southeast.  The  westerly  limit  follows  the  old  forestry  trunk  road  northwestwards  from 
Nordegg  to  the  turnoff  for  Shunda  (Baldy)  lookout.  The  trace  of  the  Brazeau  Thrust 
(Douglas,  1958)  forms  the  eastern  boundary  to  the  area. 

Highway  11,  which  crosses  the  area,  joins  the  Banff-Jasper  highway  25  miles 
to  the  west.  Eastwards,  Highway  11  leads  to  Rocky  Mountain  House  (57  miles  away) 
and  Red  Deer.  The  North-South  Forestry  Trunk  Road  crosses  Highway  11  two  miles 
west  of  Nordegg  and  connects  it  to  the  Jasper- Edmonton  Highway  to  the  north  and 
the  old  Banff-Calgary  Highway  to  the  south. 

Within  the  area,  the  Forestry  road  turning  north  from  Highway  11  at  the  Beaver 
Pools  in  Shunda  Gap  gives  access  to  the  southeast  end  of  Coliseum  Mountain.  The 
only  other  road  open  to  vehicles  is  the  trail  to  the  lookout  on  the  top  of  Shunda 
Mountain.  There  is  a  pack  trail  up  the  southwest  side  of  Coliseum  Mountain  to  the 
old  lookout. 

The  area  lies  within  the  Rocky  Mountain  Forest  Reserve.  It  is  therefore 
covered  by  4  inches  to  1  mile  aerial  photographs  and  planimetrical ly  corrected 
4  inches  to  1  mile  uncontoured  maps  of  the  tree  cover  (Forestry  Access  Maps).  A 
contoured  map  of  the  Nordegg  area  is  published  on  the  scale  of  1  inch  to  1  mile; 
1:50,000  contour  maps  of  the  surrounding  areas  are  published. 
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Stratigraphy 


Introduction 

The  stratigraphy  of  the  thesis  area  is  summarized  in  Table  1  and  in  a  geological 
map  of  the  northern  half  of  the  Brazeau  Range  (Fig.  9).  Further  detail  is  available  in 
the  publications  by  Erdman  (1950)  and  Douglas  (1956,  1958). 


TABLE  I.  Stratigraphic  succession  in  the  thesis  area. 


ERA 

PERIOD 

OR  EPOCH 

GROUP  OR 
FORMATION 

LITHOLOGY 

THICKNESS 

M 

c 

LOWER 

LUSCAR 

Sandstone, shale, coal . 

900' 

L 

s 

CRETACEOUS 

CADOMIN 

Conglomerate,  sandstone. 

25' 

o 

z 

N 

IKANASSIN 

Sandstone, shale . 

400 1 

O 

1 

JURASSIC 

FERNIE 

Conglomerate,  mudstone. 

500 1 

c 

chert,  sandstone. 

n  codm  itv 

p 

MISSISSIPPIAN 

R 

MOUNT  HEAD 

A 

U 

L 

N 

&  TURNER 

Vuggy, massive  bedded 

300 1 

A 

D 

dolomite  (undifferentiated). 

E 

L 

VALLEY 

O 

E 

Z 

SHUNDA 

Limestone,  shale,  anhydrite. 

200' 

o 

1 

dolomite  breccia . 

PEKISKO 

Limestone 

120' 

c 

BANFF 

Argillaceous  limestone, shale. 

550' 

DEVONIAN 

EXSHAW 

Shale. 

5' 

PALLISER 

Limestone,  dolomite 

800' 

ALEXO 

Sandstone,  si  Itstone,  dolomite 

200' 

:  r  torn*  » 

r  vo  i\  •  >H«*I  ft  .(?  ufcesMfl  »r .  Vori  r-  ^  hen  ••-•  •  to  <  wf 


' 


v/T  L  *•*>  •>  A.'  •  . 


- .  ..  VtIM!  CTMCJc  i 

. 


v3  J , A\ 


■. 


. 


5 


GEOMETRIC  ANALYSIS 

Determination  of  the  Orientation  of  Surfaces 

The  area  was  mapped  on  4  inch  to  1  mile  aerial  photos.  Detail  easily  correlated 
with  the  photos  made  transfer  to  the  Forestry  Access  Maps  accurate  and  simple.  Spot 
heights  have  been  determined  by  the  Topographical  Survey  on  the  peaks  of  Coliseum  and 
Shunda  Mountains  and  step  outs  from  these  were  made  by  the  writer  in  1964  by  Brunton 
compass  and  in  1965  by  altimeter. 

The  writer  concentrated  on  determining  with  a  definable  precision  the  attitude  of 
a  number  of  surfaces  throughout  the  map  area.  Practically  all  exposures  (or  stations) 
close  to  the  Alexo-Palliser,  Devonian-Mississippian  and  Banff-Pekisko  boundaries 
were  visited.  Wherever  possible  the  orientation  of  the  surface  was  established.  The 
ideal  size  and  spacing  of  the  stations  is  governed  by  the  rate  of  variation  of  the  orien- 
tation  of  the  surface.  In  areas,  such  as  fold  hinges,  where  the  surface  changes  orien¬ 
tation  rapidly,  more  stations  are  required  to  determine  the  surface  than  on  the  limbs.  The 
size  of  stations  on  fold  hinges  are  also  reduced,  for  the  statistics  used  to  analyse  the 
stations  assume  that  the  errors  in  measurement  are  randomly  distributed  about  a  mean 
direction  (with  a  certain  probability  density  function),  so  that  the  mean  orientation 
of  the  surface  should  not  change  significantly  within  the  station. 

It  is  possible  to  distinguish  three  sources  of  uncertainty  in  the  mapping  of  a 
particular  exposure  of  a  folded  surface,  uncertainty  about  the  topographic  position 
of  the  exposure,  uncertainty  about  the  stratigraphic  position  of  the  exposure  and  un¬ 
certainty  about  the  orientation  of  the  surface  within  the  exposure.  Only  the  last  of 
these  will  be  discussed  because  the  first  source  can  be  reduced  to  negligible  prop¬ 
ortions  by  sufficiently  sophisticated  survey  methods,  and  the  importance  of  the 
second  is  at  a  minimum  in  the  mapping  of  well-exposed,  well-defined  boundaries. 


eUY.lA*A  D!flT3*/0:  D 


-  .  . 


bMdv  ,9  UetsO  *  >Hq  Soji*i>  .1  ■  ■  I  01  <brf  *  «K  »•'  '  »• 

,  •-  .  .r.l.o  .v  *01*  *««.*■>':  ***-*«« 


:i  15  ,q  ,  -bo  «  ■v.loin.  •’  b  oo  I.':  .;o  :v:  v  *HV 

_ _ _ _ 


.  •  >o  it  ‘it.  ■  ■'>*'■  7<  ’  °  '  9  ‘ 

n  0,  I'lcio  *  9  ■  *  1  *'">*'*  '  "  •t9,iiiv  ”  w 

-  liio  te^norto  9  ,-W  H  .w  ^8ni(|  bid  *  *H*  •«  «<  •«  ^ 


iat^D  Gt  beiiupai  »ip  anohoU  mom  y\'  qt>»  no»+B 
,H,  <N'onoo  .  «u  .»!  °t*sf.ol  *  .oA  r  .-io  aio  *or,  «-o  mc’.-M  *o  > 


. 

■ 


..QOiq  sldieile^o  oi  bosubo,  od  noo  soiuoa  Wrt  9*  «4MM>  bo*uoiib  sd  l!5w  wodt 
•Hi  to  90001, oqm;  aifl  boo  .bodlao,  ^  b-'o  ,10-Hq.  t  XHn**3t -»u«  Xd  .00,1,0 


6 


The  importance  of  the  third  source  is  indicated  in  Mertie's  work  (1947).  Th  is  source 
does  not  include  uncertainty  as  to  whether  the  exposed  surface  is  a  sample  of  the 
covered  surface  or  has  been  affected  by  processes  such  as  hill  creep.  The  existence 
of  hill  creep  can  be  demonstrated  only  by  comparison  with  exposed  undisturbed  rocks 
and  these  can  be  used  to  determine  the  representative  orientation  of  the  surface. 

Where  hill  creep  cannot  be  demonstrated  a  decision  about  a  particular  exposure  is  a 
decision  about  a  likely  value  for  the  surface  at  this  point.  Such  decisions  belong  to  a 
dynamic  analysis  not  to  a  geometric  one. 

It  is  natural  to  assume  that  errors  involved  in  measuring  the  orientation  of  a  surface 

are  normally  distributed.  However,  if  the  data  are  to  be  treated  by  the  usual  graphical 

techniques  they  are  displayed  as  directed  lines  radiating  from  the  centre  of  a  sphere  of 

finite  radius.  Whereas  Gauss  developed  the  normal  distribution  on  a  surface  whose 

topology  could  be  ignored,  in  this  case  the  errors  involved,  often  several  degrees, 

require  the  topology  of  the  sphere  be  taken  into  account.  Fisher  (1953)  has  suggested 

how  this  might  be  done.  Suppose  the  distribution  of  possible  observations  be  on  the 

surface  of  a  unit  sphere.  He  proposed  that  the  distribution  of  elementary  errors  over 

K  c  osO 

this  surface  had  a  frequency  density  proportional  to  e  where  9  is  the  angular 

displacement  from  the  estimated  mean.  The  density  is  a  maximum  at  9  =  0  (if  K 
is  positive).  If  K  =  0  the  density  is  uniform  over  the  sphere.  If  K  is  large  and 
positive  the  distribution  is  confined  to  a  small  area  in  the  neighbourhood  of  the  max- 
imum  (i.e.  9  is  small).  Fisher's  proposal  seems  natural,  for  if  9  is  small 

cos  e  -  C 1  -  ) 

and  the  frequency  density  becomes 

-  KJ92- 
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which  is  a  2-dimensional,  normal  distribution  with  K  the  analog  of  the  inverse  of  the 

variance,  a  consequence  of  minimizing  the  effects  of  the  topology.  Integration  of 

the  function  over  the  sphere  gives  a  normalizing  constant,  c  =  - — _ 

ATTS-'ttln  K 


and  a  probability  density  function,  P  (Q), 


P(0)  = 


k 


4-7 K 


kcosB 

e 


Paleomagnetic  data  have  been  found  to  fit  this  distribution  and  Muecke  (1965) 
has  suggested,  though  not  tested,  that  multiple  readings  of  apparently  parallel  joints 
in  rocks  are  distributed  in  this  manner.  To  test  whether  this  distribution  can  be  used 
to  deal  with  replications  of  bedding  planes  from  the  same  domain,  61  dip  and  strike 
measurements  were  made  on  an  outcrop  of  Palliser  limestone  on  the  southeast  end  of 
Coliseum  Mountain. 

The  null  hypothesis  that  these  61  measurements  did  not  significantly  differ  in 
distribution  from  Fisher's  distribution  could  then  be  examined  at  the  5  per  cent  sig¬ 
nificance  level.  The  procedure  is  set  out  (for  paleomagnetic  data)  in  Watson  and 
Irving  (1957).  A  computer  calculation  (see  Appendix  4)  of  the  mean  bedding  plane 
revealed  a  dip  of  12°  to  212°  with  precision  parameter^  K  =  833.  The  poles  to  the 
bedding  planes  were  plotted  on  a  stereographic  projection  and  rotated  to  bring  the 
mean  to  the  centre  of  the  projection  (Turner  and  Weiss,  1963,  p.  67). 

Two  test  of  the  hypothesis  were  then  applied.  The  first  tests  whether  the  data  are 
significantly  divergent  from  axial  symmetry  about  the  calculated  mean  direction, 
the  second  whether  the  data  are  significantly  divergent  in  distribution  from  the 


probability  density  function,  P  (0  ) 

In  the  first  test  the  distribution  of  poles  on  the  projection  was  divided  into 
octants  and  the  significance  of  the  results  tested  against  chi-square  with  5  degrees  of 
freedom  (Lindley  &  Miller,  1953).  The  hypothesis  (that  the  data  are  a  sample 
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from  a  population  with  a  Fisher  distribution)  is  rejected  at  the  5%  level  of  significance, 
if,  on  the  assumption  that  the  hypothesis  is  true,  the  probability  is  less  than  1  in  20 
of  exceeding  the  observed  value  of  chi-square  by  chance. 
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7.6 
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7.6 
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7.6 

7.6 

sf.-$  X, 


Significance 

6.55  11.07  None 

where  fo  is  the  observed  frequency  and  fe  is  the  frequency  predicted  by  a  Fisher  dist- 

ribution . 

Watson  and  Irving  (1957)  showed  that  the  expected  frequency,  fe,  of  readings 
between  the  angular  distances  0^ ,  from  the  mean  (where  0,  is  less  than  9 )  is  for 

large  K, 

fe  =  N  (exp(-K(l-cos9j))  -  exp  (-KO-cose^)), 
where  N  is  the  total  number  of  readings.  The  hypothesis  is  again  tested  against  chi- 

X,' 


square. 

0 
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fo 
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fe 

5.7 

18.6 

16.6 

12.1 

8.0 

I 


Significance 


1.21 


3.84  None 


There  is  again  no  reason  to  reject  the  hypothesis,  so  the  distribution  of  the 
bedding  planes  in  both  tests  does  not  significantly  depart  from  Fisher's  distribution. 

It  is  worthwhile  noting  that  these  tests  can  be  more  accurately  performed  by  doing  the 
arithmetic  rather  than  by  the  graphical  method  from  Turner  and  Weiss.  The  process 

is  tedious  by  hand  but  is  suited  to  computer  techniques. 

Watson  and  Irving  (1957)  derived  several  relations  which  are  helpful  in  estimating 
in  advance  the  number  of  bedding-plane  measurements  required  at  a  particular 
station,  assuming  the  measurements  follow  Fisher's  distribution.  In  this  survey. 
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stations  were  set  up  with  one  of  two  objectives  in  mind. 

(1)  Suppose  that  it  is  required  to  determine  the  orientation  of  bedding  at  a 
particular  station  to  within  of  the  estimated  mean.  Then,  after  Watson  and 

Irving  (1957,  p.  291)  for  large  K 


-  IHC°  (  K  n)"^ 

where  N  is  the  number  of  readings,  and 


92  *s  radius  of  the  95%  cone  of  confidence. 


Fig.  2  is  a  graph  of  95  against  N  for  different  values  of  K, 


for 


K  =  1000 


K=  500 


N  =  3 


N  -  6 


N*  31,  for  o<95  =  2  1/2° 


K  =  100 

Thus,  as  long  as  K  is  known,  N  can  be  determined.  A  preliminary  analysis  showed 
that  K  estimates  usually  fell  between  200  and  700.  Ten  readings  were  thus  usually 
adequate  to  reduce  9^  to  2  1/2°,  A  fuller  discussion  is  given  in  Appendix  1 . 

(2)  It  seems  to  have  been  the  practice  (see  e.g.  Bielenstein,  1964)  in  some 
previously  published  studies  of  fractures  in  rocks  to  consider  two  groups  of  S  planes 
distinct  if  their  95%  cones  of  confidence  do  not  intersect.  Agterberg  (1964) 
pointed  out  that  the  optimal  procedure  in  two  dimensions  is  to  use  Student's  t-test. 
This  sets  up  the  hypothesis  that  both  samples  are  drawn  from  the  same  population 
and  tests  it  against  tables  of  the  distribution  of  the  variance  of  the  estimated  mean 
of  small  samples.  If  X  is  the  difference  in  the  estimated  means,  of  the  two 

samples,  with  sizes  N^,  Nj/  the  hypothesis  is  rejected  if  x 
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where  t  is  the  value  of  the  ti  distribution  (Lindley  &  Miller,  1953,  Table  3)  with 
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(N  +  N2  -  2)  degrees  of  freedom,  and  3  is  the  estimated  variance  of  the 
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the  population.  This  may  be  contrasted  with  an  analogy  of  "Bielenstein's"  procedure 


which  would  be 


X 


5, 


A 


s. '  (  - 1) 


t1  is  the  value  of  t  with  Nj  -  1  degrees  of  freedom, 


t"  is  the  value  of  t  with  ^  ■  1  degrees  of  freedom. 


In  three  dimensions,  Watson  (1956a)  has  shown  that  an  approximate  F  test  of 


the  hypothesis  can  be  derived  from  an  analysis  of  variance  analogy.  The  hypothesis  is 


rejected  if 


/V-l)  (/%'  -R  )  ">  f 


(l'J-  K,  -  KJ 


X  ,  2 
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N  =  N1  +  N2, 


R  is  the  length  of  the  resultant  vector  of  the  vectors  with  lengths  of  R^  and  R^, 
F  is  tabulated  in  Lindley  and  Miller  (1953)  Table  7. 


The  test  assumes  that  the  samples  of  sizes,  N] ,  N?,  are  drawn  from  populations  with 
the  same  K.  It  can  be  shown 
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for  two  populations  where  =  N^,  R^  =  1^'  ^  _  ^R^cos  x/2,  K  - 
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and  the  estimated  population  means  are  x  apart. 

Thus  if  K  =  500,  x  =  2°,  then  N]  12.  A  computer  program  (described  in 
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Appendix  4)  has  been  written  to  calculate  and  R.  This  is  a  sensitive  method 

of  measuring  small  angular  unconformities.  The  technique  is  presented  in  Appendix  3. 

To  summarize,  no  more  than  10  readings  are  usually  required  at  a  station  to 
determine  the  orientation  of  a  surface  at  that  station  to  within  2  1/2  degrees,  and 
about  12  are  required  to  distinguish  it  from  another  station  whose  mean  differs  by 
more  than  2  degrees. 


Determination  of  Fold  Geometry 

Folds  are  usually  described  by  specifying  their  axes  and  the  orientation  of 
their  axial  planes.  This  thesis,  however,  presents  the  data  as  structure  contour  maps 
on  parts  of  three  surfaces,  the  base  of  the  Pa  Miser,  the  base  of  the  Banff  and  the  base 
of  the  Pekisko.  It  will  be  shown  that  to  construct  the  contour  diagram  the  folded 
surface  must  be  divided  into  cylindrical  portions  and  the  axes  of  these  portions  calcu!-= 
ated.  From  these  axes,  the  axial  traces  and  hence  the  axial  plane  of  the  main  fold  can 
be  constructed.  There  is  no  loss  of  information  involved  in  this  mode  of  presentation. 
In  fact,  the  contour  diagram  presents  considerably  more  information  than  the  convent¬ 
ional  method  of  description .  For  instance  it  contains  sufficient  information  to  enable 
the  drawing  of  cross  sections  of  the  folded  surface  in  any  direction. 

Construction  of  structure  contour  maps 

The  contour  maps  on  the  three  surfaces  have  been  constructed  almost  indep¬ 
endently  of  one  another;  very  little  of  the  information  used  to  construct  one  map  has 
been  used  in  the  construction  of  any  of  the  other  maps.  This  restriction  causes  more 
extrapolation  into  covered  areas  or  areas  of  poor  exposure  but  it  avoids  the  effects  of 
the  restrictive  assumptions  made  in  the  downward  and  upward  projection  of  information. 
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Oertel  (1962)  has  given  a  qualitative  demonstration  of  the  importance  of 
symmetry  axes  in  extrapolating  data.  It  is  useful  to  demonstrate  this  technique  more 
formally.  Suppose  part  of  a  fold  can  be  shown  to  be  cylindrical  (and  it  seems  a  gen¬ 
eral  assumption  that  any  fold  may  be  divided  into  cylindrical  portions,  see,  for 
example.  Turner  and  Weiss,  (1963,  p.  163))  Consider  a  small  part  of  the  generatrix 
(or  fold  axis)  of  that  surface  at  a  height,  P.  The  locus  of  all  points,  P  ,  lying  on  the 
contour  line  x  feet  below  P  is  such  that  the  horizontal  distance  separating  P  from  P^ 
parallel  to  the  trend  of  the  fold  axis  is  Wof.X  where  the  plunge  of  the  axis  is  . 
This  construction  has  the  advantage  that  it  concentrates  uncertainty  about  the 
form  of  the  surface  in  uncertainty  about  the  position  of  the  contour  line,  the  end 
result  of  the  construction  (Fig.  3) 

Calculation  of  the  Generatrix 

Fleuty  in  a  review  paper  (1964)  pointed  out  that  "the  concept  of  the  calculated 
fold  axis  is  fundamental  to  much  modern  structural  analysis".  Techniques  of  calcul- 
ating  the  axes  of  cylindrical  folds  have  received  a  good  deal  of  attention  but  no 
recent  additions  have  been  made  to  the  two  common  graphical  techniques,  using  pi 
and  beta  diagrams,  for  nearly  forty  years.  However  dissatisfaction  with  these  tech¬ 
niques,  particularly  the  beta  intersection  technique,  had  been  growing. 

Ramsay  (1964)  pointed  out  that  beta  diagrams  have  the  disadvantages  that  the 
number  of  beta  axis  determinations  increases  very  rapidly  with  the  number  of  obser¬ 
vations  rntil  graphically  the  data  becomes  unwieldy;  errors  in  data  lead  to  spurious 
beta  axes;  and  that  the  diagrams  do  not  yield  as  much  information  about  the  shape  of 

folds  as  pi  diagrams. 

Muecke  (1965)  has  shown  how  the  first  disadvantage  can  be  overcome  using  a 
computer  program.  If  0  ] ,  rry  nry  n^  are  the  direction  cosines  of  the 
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outward  directed  normals  to  2  bedding  planes  then 


(I3/  ni^,  n^)  are  the  direction  cosines  of  the  line  of  intersection. 

To  find  the  mean,  he  assumed  that  the  distribution  of  the  beta  intersection 
determinations  is  a  Fisher  distribution  round  the  mean  direction.  He  also  used  this 
distribution  to  set  up  a  test  of  whether  the  fold  is  cylindrical,  using  the  precision 
parameter  and  implicitly  suggesting  the  use  of  Watson's  tables  for  random  K  (1956b). 
Unfortunately  this  elegant  solution  fails  since  it  can  be  shown  that  the  intersection 
does  not  have  a  Fisher  distribution  (Appendix  2),  The  actual  form  of  the  distribution 
of  the  intersections  seems  extremely  complex  and  it  is  not  likely  that  further  work 
will  give  any  convenient  statistical  tests  (MacGregoijl  963).  Ramsays  second  and 
third  objections  still  stand  however,  and  impair,  the  usefulness  of  the  beta  inter- 
section  technique. 

The  pi  diagram  also  has  its  disadvantages;  there  is  no  clear  criterion  for  deciding 
the  "best  fit"  axisjdrawing  the  best  fit  great  circle  by  eye,  as  is  usual,  leads  to  incon¬ 
sistency  and  error,  for  equal  distances  on  the  projection  do  not  represent  equal 
angular  distances  on  the  sphere;  there  are  finite  errors  in  plotting  readings  and  con-* 
structing  the  diagrams  -  the  maximum  likely  error  in  plotting  a  reading  is  probably 

about  2  degrees. 
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The  contouring  of  pi  diagrams  is  discussed  at  some  length  in  Appendix  5  but 


even  supposing  the  problems  associated  with  contouring  pi  diagrams  had  been  solved 
by  the  method  proposed  there,  two  disadvantages  of  the  method  remain.  Both  these 
disadvantages  and  the  contouring  problem  can  be  overcome  if  the  generatrix  is  cal¬ 
culated  mathematically.  It  is  also  desirable  to  test  whether  the  fold  is  cylindrical. 

If  the  fold  is  not  cylindrical  the  direction  of  the  axis  calculated  as  described  below 
is  meaningless.  To  perform  this  test,  the  effects  due  to  uncertainty  about  the  position 
of  the  pole  must  be  separated  from  those  due  to  uncertainty  about  the  fit  of  the  poles 
themselves  to  a  cylindrical  surface.  This  requires  a  prior  estimate  of  the  variance  of 
the  individual  readings,  obtained  by  taking  sufficient  bedding  plane  readings  from 
the  population  of  bedding  plane  readings  at  each  station  to  give  an  adequate  estimate 
of  the  variance  of  the  population. 

Suppose  then  that  p  stations  (or  outcrops)  on  the  fold  have  been  sampled  and 
th 

that  the  i  of  these  p  samples  contains  N.  readings  and  has  a  resultant  with  direction 
cosines  (I.,  m.m  n.)  of  length  R..  Then,  providing  that  the  samples  are  drawn  from 

iii  i 

populations  having  Fisher  distributions,  following  Watson  (1960)  it  can  be  shown  by 
an  analysis  of  variance  analogy  that  the  estimate  of  the  direction  of  the  fold  axis  is 
the  eigen  vector  corresponding  to  the  least  eigen  value  of  the  matrix  V, 

I."  2X'I>,’ 


n  •  m 


R;f\;  nv,  n 


providing  that  the  precision  parameters  of  the  p  populations  are  the  same  and  subject 
to  a  number  of  approximations.  These  assume  that  K  is  large  and  that  all  the  N.  are 
large.  The  effect  of  departures  from  these  assumptions  is  difficult  to  unravel  and  Watson 

does  not  discuss  this. 


•A  il  r4tgn*-l  emo2  id  *t  ?mmgoib  iq^o  gnhuo+ndo  VT 


o  si  •?  r;  dfji  $  srlt  ti  s  >odi9vo  noo  rn>*  owq  901  -'-..we  It  b«o  '  ^°2'b 


-  q  ri  uorl.  '  ci  i  b  .:•■»•>  *  «.<1  r*»  «  *'  .r-Wfq-oT  in  m 

.5 


D  "  -••••>  '  "  "  '!fr 


,  roitoiuqo  H  k  >UDV  9rit  *0 


,  ■•■.:,  •  ii\y  ADtlo  3.  .  oH  300,  ■;  btfcn  >1  ?ni  .*00  zeqroi*  <:  w  ’ 


;  V  .  :  -  !  •'  I  vt  .*  ib  Mi  »0  910  -.12-  >fi  •  OI0.1C  >'-  IDi  DV  O  .-O,.  OOD  - 


■  'v  ;0  4  :  t  ■*■■'•■  •  '  ■  •  ‘  — *  '  *'}  ■'  ' 


' 


i  ’■»?  t  .  J  9fh  -  -.o  enoi Jo,iuqoq  q  ei  »*1o  nst&'nov>q  ncia'oarq  -?rl1  for!  g.ubWoiq 


Ho  *ofh  bno  so^ol  X  loHl  ‘imuWD  se^rl",  .zftoitomixoiqqD  <0  lodmun  <  o 
idvo^nu  0‘  |iurO:ib  •  lo'tQmu..  *29  Hi  no;)  .-.yiut'oq^b  V  9  '.:  >t 


.ii  ctuviib  1o<i  ,*  »ob 


15 


A  test  of  coplanarity  is  given  by 


^  ,  •  (least  eigen  value  of  V) 


T 


—  2- 


The  principal  radii  9^,  the  cone  of  confidence  about  the  axis  are  given  by 


where  A  ,  A.  ,  A  .  are  the  maximum,  intermediate  and  minimum  latent  roots 
max  int  mm 

of  V.  Watson  (1960)  goes  on  to  show  that  for  unequal  K,  W  should  be  computed. 
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The  test  of  coplanarity  is 
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The  direction 
eigen  value. 


of  the  normal  to  the  best  fit  circle  is  the  eigen  vector  of  the  least 
Confidence  limits  about  this  direction  are  given  by 
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where  the  A  are  the  corresponding  eigen  values  of  W. 

Notice  that  in  W  the  station  averages  are  weighted  not  only  by  their  length 
but  by  their  precision  parameters.  (The  best  estimate  of  K.  is 
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Outstanding  problems  with  the  general  use  of  W  to  calculate  B  axes  lie  in  the 
unexplored  effects  of  Watson's  approximations  and  the  unverifiable  character  of  the 
results.  The  calculations  are  to  problem;  an  IBM  7040  computer  can  be  persuaded  to 
compute  an  average  size  W  in  about  half  a  minute. 

How  large  should  N.  be  before  the  approximations  Watson  uses  are  valid? 
Watson  (1960)  considers  an  example  with  N.  as  low  as  9.  In  a  later  paper  (1965) 
he  treated  data  where  N.  has  the  value  1  and  K  for  the  populations  is  known.  If 
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N.  is  1,  V  takes  the  form,  U 
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This  is  the  matrix  Scheidegger  (1965)  derived  by  a  different  method  using  a  different 
statistical  distribution.  He  used  it  to  solve  an  analogous  problem-- "the  determin¬ 
ation  of  the  'best'  axis  from  each  of  a  series  of  axes  (the  normals  to  the  bedding  planes 

or  the  long  axes  of  sand  grains)". 

In  Watson's  analysis,  the  test  for  coplanarity  is 
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The  results  of  this  test,  he  commented,  have  no  sound  probability  basis  but  are  use¬ 
ful  in  describing  the  distribution.  Somewhat  arbitrarily  the  least  permissible  value 
of  N.  in  this  study  has  been  set  at  10. 

How  may  the  fit  of  great  circle  be  tested?  Creer  (1962)  argued  that  the  fit  of 
the  best  great  circle  is  satisfactory  (at  the  5%  level?)  if  the  great  circle  passes  through 
the  5%  error  circles  around  the  mean  directions  of  the  p  populations.  This  is  not 
however  a  sufficient  condition.  The  probability  that  the  true  means  at  p  stations 

simultaneously  lie  inside  the  5%  error  circles  around  their  estimated  means  is 

P 

(0.95)  (less  than  0.95).  So  that  the  radii  of  the  circles  to  be  used  in  the  check 
is  properly  smaller  than  the  5%  circles.  But,  since  the  distribution  of  probability 
around  the  true  mean  varies  with  the  precision  parameters  of  the  population,  this 
whole  approach  is  a  little  oversimplified.  The  eventual  solution  of  the  problem  of 
finding  an  independent  check  may  lie  in  methods  of  the  type  suggested  by  Selby 
(1965). 

Discussion  of  procedure 

Because  Watson's  method  gives  an  objective  test  of  the  domain  of  a  cylind¬ 
rical  fold,  the  matrix,  W,  has  been  used  to  calculate  the  generatrix  of  cylindrically 
folded  portions.  A  computer  program  to  perform  these  calculations  is  given  in 

Appendix  4. 

After  three  neighbouring  stations  which  lie  on  the  same  cylindrical  fold  have 
been  found,  neighbouring  stations  were  added  to  the  domain  until  the  addition  of 
any  further  station  gave  insignificant  results.  In  this  way  nearly  all  stations  within 

the  thesis  area  were  assigned  to  a  domain. 

Some  stations  fit  in  two  or  more  adjacent  domains,  but  in  general,  surfaces 
have  been  divided  into  discrete  domains  by  the  analysis.  The  positioning  of  the 
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boundaries  of  domains  is  rather  arbitrary.  Obviously  they  lie  somewhere  between  the 
last  station  within  the  domain  and  the  first  outside  it,  but  where  is  difficult  to  say. 
Conventionally  they  were  placed  at  the  halfway  mark. 

It  is  possible  that  dividing  the  surfaces  into  domains  is  an  artificial  process 
much  like  the  statistician's  display  of  a  continuous  variate  as  a  histogram--a  comput¬ 
ational  convenience.  It  may  have  the  effect  of  concentrating  changes  in  shape  of 
the  fold  along  the  margins  of  the  domains  in  an  unrealistic  manner.  Although  the 
analysis  reveals  a  fold  axis  single-valued  throughout  the  domain,  it  is  possible  that 
the  axis  varies  systematically  but  insignificantly  across  the  domain.  The  analysis 
takes  no  account  of  the  relative  positions  of  the  stations  in  space.  This  effect 
might  be  examined  by  dividing  large  domains  into  a  number  of  subdomains  and  calcu¬ 
lating  the  axes. 

Since  the  analysis  takes  no  account  of  the  relative  positions  of  the  stations 
these  might  be  expected  to  provide  some  additional  information  about  the  shape  of 
the  fold.  This  information  can  be  used  in  two  ways.  (1)  If  the  axis  is  designated 
by  the  most  probable  value  given  by  the  analysis,  then  the  information  can  be  used  to 
construct  the  contours  and  thus  the  shape  of  part  of  the  fold.  (2)  There  is  no  reason 
to  assume,  however,  that  when  the  new  information  is  considered  the  probable  value 
of  the  fold  axis  is  still  the  same.  It  might  be  supposed,  for  instance,  that  a  value  of 
the  fold  axis  should  be  chosen  that  keeps  the  contour  lines  as  smooth  as  possible  as 
this  is  the  most  "natural"  shape  for  a  contour  line.  In  practice  the  procedure 
adopted  compromises  between  these  two  extremes.  Contours  have  been  constructed 
by  displaying  the  limits  the  calculated  fold  axis  places  on  them  and  then  drawing 
contour  lines  to  pass  within  these  limits  subject  to  the  condition  that  the  fold  axis 
is  single  valued . 
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In  detail,  the  program  outputs  a  most  likely  value  for  the  fold  axis  together 
with  the  lengths  and  directions  of  the  semi-axes  of  the  elliptical  cone  of  confidence 
about  the  axis.  Then  it  has  been  shown  that  the  next  lowest  contour  line  is  xcota 
feet  away  (  p.  )  down  the  trend  of  the  fold  axis.  Nearly  all  the  fold  axes  en¬ 
countered  dip  gently.  This  has  the  effect  of  producing  an  ellipse  of  confidence 
whose  projection  down  the  fold  axis  has  negligible  thickness.  It  can  be  approximated 
by  a  straight  line.  The  contour  line  must  then  pass  between  two  points  x  cot(a-a'), 
x  cot(a+a')  from  the  point  of  projection  down  the  trend  of  the  fold  axis  where  a1 
is  the  angle  the  major  semi  axis  of  the  ellipse  of  confidence  subtends  at  the  apex  of 
the  cone.  If  the  value  of  x  is  certain  (see  p.  )  then  the  only  effect  of  x  is  to 
increase  the  area  of  uncertainty  as  its  value  increases,  a  natural  consequence  of 
increasing  distance  from  the  point  of  projection. 

Suppose  the  compromise  leads  to  a  value  a"  for  the  plunge  of  the  fold  axis. 

Then  the  contour  line  250  feet  below  the  constructed  contour  line  is  (250  cota"  + 
x  cota"  feet),  (250  cota"  +  x  cot(a-a')),  (250  cota"  +  x  cot(a+a'))  from  the  point  of 
projection.  The  limits  on  the  lines  are  displayed  on  the  contour  diagram  of  the  top 
of  the  Alexo  formation.  On  the  other  diagrams  these  limits  have  been  omitted  for 
drafting  convenience. 

As  there  is  little  exposure  north  east  of  the  trace  of  the  highest  splay  of  the 
Brazeau  Range  fault  there  was  little  point  in  extending  contours  across  it.  On  the 
map  of  the  top  of  the  Palliser  this  has  made  the  contouring  on  the  Palliser  below  the 
thrust  plane  easier  to  appreciate.  There  is  of  course  no  discontinuity  in  the  contours 
at  the  trace  of  a  fault  plane.  The  discontinuity  exists  only  where  the  contour  line 
and  the  trace  of  the  fault  plane  are  at  the  same  height. 
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KINEMATIC  ANALYSIS 
The  Northern  Part  of  the  Brazeau  Range 

To  make  a  reasonable  kinematic  analysis,  the  rest  of  the  major  structure  of 
which  the  thesis  area  is  part  must  be  considered.  Material  for  this  has  been  taken  from 
Douglas  (1956,  1958)  and  Erdman  (1950). 

The  thesis  area  lies  within  the  Brazeau  thrust  sheet,  towards  the  northwest 
end  of  the  Brazeau  Range.  The  Lower  Mesozoic  and  the  Palaeozoic  strata  that  make 
up  the  Range  have  been  brought  to  the  surface  along  the  Brazeau  Range  fault,  a  splay 
of  the  Brazeau  thrust.  These  strata  are  folded  into  the  Brazeau  Range  anticline.  On 
the  Chungo  Creek  sheet  (Douglas,  1958)  the  anticline  plunges  northwestwards  and 
disappears  ten  miles  northwest  of  Shunda  Gap.  A  southeast  limit  to  the  anticline  can 
be  seen  near  the  head  of  Lundine  Creek  on  the  margins  of  the  Alexo  and  Saunders 
map  areas.  "Further  to  the  southeast  the  Brazeau  Range  anticline  plunges  below  the 
surface  as  a  wide  plate  wi  th  no  definite  axial  trace."  (Erdman,  1950,  p.  23).  The 
total  length  of  the  anticline  is  twenty-four  miles. 

Befween  Shunda  Creek  and  the  North  Saskatchewan  River  the  strike  of  the 
axial  surface  of  the  anticline  is  approximately  135  degrees.  North  of  Shunda  Gap 
the  strike  changes  gradually.  Where  the  anticline  dies  out  on  the  Chungo  Creek  sheet 
the  strike  is  about  100  degrees. 

The  Brazeau  Range  anticline  is  a  multiple-plunging  structure.  There  are 
two  major  culminations  north  of  the  North  Saskatchewan  River,  one  three  miles  northwest 
of  the  river  where  the  top  of  the  Palliser  is  at  a  height  of  7,500  feet,  and  the  other 
the  culmination  on  Coliseum  Mountain  mapped  in  this  thesis,  where  the  top  of  the 

Pa II iser  is  at  6,100  feet . 
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The  strata  lying  between  the  Brazeau  Range  fault  and  the  Brazeau  thrust 
itself  are  folded  into  the  Swale  Creek  syncline.  North  of  the  more  southerly  culmin¬ 
ation  in  the  Brazeau  Range  anticline,  strata  on  the  northeast  limb  of  the  anticline  rest 
on  successively  older  strata  in  the  southwest  limb  of  the  syncline:  Luscar  strata,  for 
instance,  rest  on  Brazeau  just  north  of  the  culmination,  and  on  Blackstone  where  the 
Range  dies  out  in  the  extreme  northwest. 

A  kinematic  analysis  of  the  northern  half  of  the  Brazeau  Range  should 
explain  the  general  northwesterly  plunge  of  the  Brazeau  Range  anticline  and  the  arcuate 
strike  of  its  axial  surface.  Fig.  10  is  a  simplified  model  explaining  how  the  plunge 
of  the  anticline  and  the  arcuate  strike  follow  from  the  initial  assumption  that  there  is  a 
decrease  in  transport  northwestwards  along  the  Brazeau  Range  fault.  The  interpret¬ 
ation  of  folding  in  the  model  (after  Gwinn,  1964)  has  already  been  put  forward  by 
Douglas  (1950)  and  Hume  (1957).  Though  details  of  the  theory  have  been  criticized 
by  Scott  (1954)  and  Fox  (1959)  their  criticisms  were  not  aimed  at  the  "thin-skinned" 
interpretation  of  Foothills  tectonics  so  much  as  at  Douglas'  choice  of  several,  minor 
zones  of  decollement.  While  there  is  still  some  dispute  over  possible  decollement 
zones  within  the  sedimentary  column,  the  uncoupling  of  the  basement  from  the  over- 
lying  deformation  is  confirmed  by  seismic  evidence  and  there  now  seems  little  doubt, 
that  at  (or  just  above)  the  contact  of  the  PreCambrian  basement  complex  with  the 
overlying  sediments  there  is  a  major  zone  of  decollement  under  much  of  the  Foothills 
and  Eastern  Ranges  of  the  Rocky  Mountains. 
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The  Thesis  Area 

One  and  a  half  miles  southeast  of  Shunda  Gap,  the  Pa lllser  is  exposed  in 
the  valleys  of  the  streams  that  flow  northeast  to  join  Shunda  Creek.  The  contour  map 
on  the  Palliser  top  (Fig.  7)  shows  it  is  folded  into  two  anticlines— the  equivalent  of 
the  Brazeau  Range  anticline  and  a  smaller  anticline  northeast  of  this.  The  intervening 
syncline  is  almost  isoclinal  towards  its  apex,  the  limbs  dipping  towards  one  another  at 
80  degrees  in  the  extreme  southeast  of  the  thesis  area  (Fig.  7). 

It  is  possible  to  trace  Palliser  continuously  through  the  same  two  anticlines 
and  tight  syncline  on  the  hillside  looking  into  Shunda  Gap  from  the  southeast.  The 
syncline  in  the  Palliser  is  isoclinal,  both  limbs  dipping  to  the  southwest  at  about  65 
degrees.  In  the  overlying  Banff  strata,  the  syncline  is  less  tight  with  uninverted  dips 
from  30  to  40  degrees  in  the  limbs,  evidence  of  structural  disharmony. 

The  next  exposures  to  the  northeast  are  in  the  cuts  along  the  railway  line 
in  Shunda  Gap.  The  syncline  and  anticline  northeast  of  the  Brazeau  Range  anticline 
are  no  longer  clearly  identifiable.  The  section  shows  a  series  of  fault  slices  of  over- 
turned  strata,  the  top  of  the  Palliser  and  most  of  the  Banff  Formation  are  thrust  in 
contact  with  Mount  Head  and  Fernie  rocks  which  themselves  lie  in  thrust  contact 
with  Cadomin  and  Luscar. 

A  thrust  has  been  inserted  between  the  northeast  dipping  Palliser  on  the 
northeast  limb  of  the  Brazeau  Range  anticline  and  the  overturned  southwest  dipping 
Palliser,  even  though  there  is  no  space  problem  on  the  line  of  section.  The  map  of 
the  Palliser  top  (Fig.  7)  shows  this  fault  is  necessary  bgth  to  the  southeast  and  north- 

west  of  the  Gap. 

Douglas'  (1956)  Section  CD  which  here  lies  almost  along  the  railway 
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line,  inserted  an  anticlinal  hinge  below  the  fault,  implying  that  the  fault  had  locally 
migrated  to  the  back  limb  of  the  Brazeau  Range  anticline.  A  look  at  the  contour  map 
on  the  top  of  the  Alexo  (Fig.  6)  shows  that  this  may  have  oversimplified  the  problem. 
On  the  southeast  side  of  Shunda  Gap  the  fold  plunges  in  a  direction  slightly  north  of 
northwest,  on  the  northeast  side  the  crest  of  the  fold  in  the  Alexo  is  not  exposed  but 
appears  offset  considerably  to  the  northeast  with  the  fold  axis  dipping  into  the  Gap. 
Strike  lines  drawn  on  the  fault  trace  indicate  that  the  fault  plane  has  probably  changed 
attitude  rapidly  across  the  Gap  from  a  gentle  dip  to  slightly  west  of  south  to  a  steeper 
dip  to  slightly  south  of  west.  This  sharp  change  shows  up  in  the  domains  constructed 
on  the  contour  map:  southeast  of  the  Gap  the  domains  share  two  stations,  and  the 
change  of  dip  and  strike  is  gradual,  whereas  these  domains  share  no  stations  with 
the  domain  to  the  north  of  the  Gap.  The  minor  fault  appears  to  displace  the  Brazeau 
Range  anticline  and  is  probably  later  than  the  initiation  of  the  fold. 

No  further  discussion  of  the  railway  section  is  necessary  till  the  Luscar 
exposure  is  reached.  This  seems  to  be  overturned  from  the  altitude  of  2  examples  of 
truncated  cross  bedding  and  it  is  probable  that  the  Cadomin  conglomerate  above  it 
belongs  to  the  same  thrust  slice.  Careful  measurement  shows  that  there  is  room  in  the 
thrust  slice  above  the  Cadomin  only  for  most  of  the  Fernie,  so  the  Nikanassin  is  cut 
out  of  the  section.  The  Cadomin  and  Luscar  slice  is  then  part  of  the  overturned  north¬ 
east  limb  of  the  minor  anticline.  The  Palliser  and  Banff  slice  comes  from  the  south¬ 
west  limb  (also  overturned)  of  this  anticline  and  there  is  sufficient  room  for  the  syn¬ 
cline,  almost  plunged  out,  above  the  fault.  On  the  northwest  side  of  the  Gap  the 

syncline  has  disappeared . 


. 

-  arh  b  rrsmav©  ©vori  vom  fcifto'  icrll  ewoHa  (6  .  qH)  OxelA  9^1  vo  qot  ©rit  no 


i  i  bs  :■  q  'i  ox‘i  -  a  -1  Ir*  9rH  ’o  ,  v:o-  r  •  ©hi?  von-nor  *-v:i  >ir  i 


.  !  )  '  :■  2 

•iorta  yfcbdoiq  ^  <«  i!  ji  f  ©Hi  Trrit  ©t&aibnt  sooit  lluD^  m" 

' 

. 


, 


.  ••  vi  r  '  r  .  ■  ,  i  ;tacfc  i'q  ?!  be  ■:•  ?.m  l:>Hn  ?  ognoSl 


■  , 


©fit  n  <  t  ••;  e  srH  1  *ft;  ;  yoHe  train©;  •;>  ;  i  '  ;  d  )  *;  m  -r  .  c  t  *<J 

i 

i  ,  ■  II -  i  o?  \  ■  >1  ,  i rro :3c  .  r...t  Sii?  ti  nrt 


■ 

ri:t  '  »  i  rlt  moH  .  rioo  9  •>!?  •  nofi  bno  iwiHflfl  ftriT  .on  ioitno  lortim  ©Ht  I  nr© 

j 


24 


Douglas  mapped  the  highest  minor  thrust  as  dying  out  in  the  Palliser  a  mile 
northwest  of  the  railway  line.  In  this  thesis  it  is  interpreted  as  running  between  the 
overturned  and  upright  strata  on  the  northeast  flank  of  Coliseum  Mountain  till  it 
dies  out  under  the  minor  culmination  shown  on  the  contour  map  of  the  top  of  the  Pall¬ 
iser. 

Observation  suggests  that  there  is  no  room  problem  in  the  Banff  strata 
where  Douglas  maps  a  lower  thrust  fault  nor  are  there  signs  of  the  trace  of  such  a 
fault.  But  along  the  northeast  side  of  Coliseum  Mountain  the  strata  between  the  top 
of  the  Pekisko  and  the  base  of  the  Fernie  are  almost  invariably  disturbed  and  in  a 
number  of  sections  there  is  a  definite  room  problem.  This  lower  thrust  fault  moves 
Palliser  and  Banff  strata  onto  Mount  Head  and  Fernie  strata  in  the  railway  section  and 
probably  terminates  against  the  highest  minor  fault. 

Both  minor  faulting  and  minor  folding  within  the  thesis  area  seem  to  give 

similar  results  kinematically,  they  transfer  some  of  the  horizontal  portion  of  the  total 

\ 

transport  along  the  Brazeau  Range  fault  to  the  vertical  portion.  This  is  more  easily 
seen  in  the  case  of  the  faults.  The  splaying  of  the  Brazeau  Range  fault  distributes 
total  movement  along  it  over  a  number  of  minor  faults.  As  these  lie  above  the  main 
thrust  they  have  a  steeper  dip  to  the  west  than  the  main  thrust  and  a  greater  proportion 
of  the  total  movement  on  them  is  vertical  (assuming  that  the  movement  on  these 
faults  is  dominantly  dip-slip). 

The  argument  in  the  case  of  the  minor  folding  is  less  clear  since  the 
alternative  style  of  deformation  is  not  as  exactly  specified.  However  by  comparison 
with  the  section  to  the  southeast  where  there  are  no  minor  structures  on  the  northeast 
limb  of  the  main  anticline,  we  can  argue  that  the  effect  of  the  minor  fold  is  to  reduce 
the  "apparent  shortening"  (or  horizontal  transport)  on  the  Brazeau  Range  Fault.  We 
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may  also  note  the  sequence  northwestwards  along  strike  of  no  minor  structures, 
minor  folding  and  then  splaying  of  the  Brazeau  Range  fault  and  argue  that  the  minor 
faulting  is  a  more  intense  expression  of  the  same  kind  of  deformation  that  caused 
minor  folding. 

The  result  of  the  minor  faulting  and  folding  is  the  local  reversal  of  the 
plunge  of  the  Brazeau  Range  anticline.  The  contour  map  on  the  top  of  the  Palliser 
shows  that  this  surface  is  carried  up  on  the  minor  faults  to  reach  a  culmination  slightly 
southeast  of  the  peak  of  Coliseum  Mountain  before  the  dying  out  of  two  of  the  splays 
of  the  fault  cause  it  to  resume  its  downward  plunge.  The  contour  map  on  the  top 
of  the  Banff  shows  a  culmination  further  to  the  northwest.  The  northwesterly  displace¬ 
ment  of  the  culmination  is  probably  due  to  the  smaller  curvature  of  the  Banff  surface — 
it  is  a  greater  distance  from  the  fault  plane  than  the  Palliser  surface. 

The  change  of  fold  profile  along  strike  is  also  worth  noting.  The  top  of 
the  Banff  is  almost  vertical  in  the  northeast  limb  of  the  Brazeau  Range  anticline 
above  the  highest  splay  of  the  Brazeau  Range  fault.  When  this  dies  out  the  north¬ 
east  limb  of  the  anticline  slackens  dip  and  the  hinge  of  the  fold  migrates  five  hundred 
feet  down  the  northeast  limb  from  near  the  crest.  This  should  be  explained. 

It  has  already  been  argued  that  the  initiation  of  folding  was  prior  to  the 
initiation  of  minor  faulting.  The  character  of  the  folding  might  be  expected  to  local¬ 
ize  the  minor  faulting.  This  argument  has  already  been  used  to  show  that  minor  faulting 
evolved  from  minor  folding.  Unfortunately  the  thrust  splays  show  no  signs  of  dis¬ 
appearing  into  minor  folds  at  their  northwesterly  ends  and  to  preserve  the  hypothesis  a 
"conventionalist  twist"  (Popper,  1957)  has  to  be  introduced.  Once  minor  faulting 
has  taken  place  it  will  be  assumed  that  it  spreads  faster  than  minor  folding  and  event¬ 
ually  moves  into  areas  where  minor  folding  has  yet  to  be  initiated  and  so  relieves 
the  necessity  for  minor  folding.  As  the  major  fold  has  been  presumed  in  the  kine- 
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mafic  analysis  to  have  been  produced  by  "up  and  over"  movement  along  the  Brazeau 
Range  fault,  it  is  likely  that  movement  along  a  minor  fault  will  cause  some  modifi¬ 
cation  by  'drag"  of  the  profile  of  the  major  fold.  Drag  causes  a  hinge  to  develop 
above  the  fault  plane  and  explains  the  migration  of  the  hinge  down  the  northeast 
limb  to  just  above  the  Coliseum  Fault.  There  remains  the  problem  of  why  minor 
faulting  and  folding  developed  where  they  did.  There  are  at  least  two  plausible 
solutions  to  this  problem. 

The  geological  map  of  the  northern  half  of  the  Brazeau  Range  shows  that  the 
northward  decrease  in  transport  on  the  Brazeau  Range  fault  is  apparently  reversed  in 
Shunda  Gap.  On  both  sides  of  the  Gap  the  Luscar  strata  lie  on  Brazeau  strata,  in 
the  Gap  itself  they  lie  on  Wapiabi  strata.  Douglas'  (1958)  map  suggests  that  there  is 
a  local  doming  of  the  axis  of  the  (Swale  Creek)  syncline.  Such  a  doming  would  force 
doming  of  the  structures  above  and  behind  it  and  cause  some  horizontal  movement  to 
be  transferred  to  vertical  movement. 

Alternatively,  the  sharp  change  in  strike  of  the  fold  across  the  Gap  might 
be  taken  as  evidence  that  there  is  a  sudden  increase  in  resistance  to  horizontal  move¬ 
ment  along  the  Brazeau  Range  fault  and  this  causes  an  increase  in  the  vertical  portion 
of  total  transport.  Such  an  increase  might  be  due  to  changes  in  the  properties  of 
the  horizon  of  decollement.  The  section  drawn  by  Shaw  (1963)  (which  presumably 
makes  use  of  Imperial  Oil's  unpublished  geophysical  data)  shows  that  the  Brazeau 
Range  is  the  most  northeasterly  thrust  slice  in  the  Disturbed  Belt  of  Central  Alberta 
with  a  horizon  of  decollement  below  the  top  of  the  Cambrian.  The  Devonian  and 
Cambrian  strata  under  the  Stolberg  anticline,  for  instance,  to  the  northeast  are 
believed  to  be  undisturbed.  This  change  of  horizon  of  decollement  might  be  due  to  a 
change  in  the  properties  of  the  horizon  in  which  this  decollement  takes  place  — 
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perhaps  a  facies  change  from  shale  to  limestone.  This  kind  of  facies  change  on  a  larger 
scale  might  reasonably  account  for  the  northwestward  plunge  of  the  whole  Range. 

In  summary  the  assumptions  the  kinematic  analysis  makes  are  listed  as 
they  were  developed. 

1)  The  analysis  requires  a  "thin-skinned"  model  for  the  tectonics  of 
the  Disturbed  Belt. 

2)  Movement  on  the  Brazeau  Range  fault  dies  away  northwestwards  with 
the  Brazeau  Range  anticline.  The  decrease  in  movement  may  be  due  to  a  change 

in  the  mechanical  properties  of  the  layer  of  decollement. 

3)  The  sequence  of  events  within  the  area  was  probably  --  major  faulting 
and  folding  initiated  more  or  less  contemporaneously  (the  folding  is  caused  by  movement 
up  the  fault).  Then  minor  folding  'followed  closely  by  minor  faulting1.  The  formation 
of  minor  folds  probably  ended  before  any  of  the  other  forms  of  deformation,  though 

the  major  folds  continued  to  be  modified  by  movement  on  minor  and  major  faults. 

The  dynamic  analysis  begins  the  testing  of  these  assumptions.  They  imply 
much  about  the  mechanical  properties  of  materials  which  behaved  in  this  way.  A 
solution  to  the  dynamic  analysis  would  consist  of  finding  a  material  with  properties 
consistent  with  present  knowledge  of  the  behaviour  of  rocks  under  stress  and  finding 
a  stress  field  history  consistent  with  present  knowledge  of  a  likely  stress  history  which 
would  cause  such  movements  as  have  been  postulated.  This  would  be  a  fairly  rigorous 
test  of  the  kinematic  analysis. 
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DYNAMIC  ANALYSIS 

Shaw  (1963)  has  compiled  a  useful  bibliography  of  the  rather  scanty 
literature  on  Foothills  tectonics  up  to  1961 .  It  is  enough  to  show  that  geologists 
are  occupied  with  a  model  of  a  rock  mass  slipping  towards  the  northeast  over  a  rough 
basement.  This  kind  of  model  immediately  suggests  experiments  with  a  sand  box,  a  box 
f  i  1 1  ed  with  layers  of  sand,  one  end  of  which  is  a  movable  plunger  which  is  advanced 
to  force  the  material  northeastwards. 

Link  (1949)  uses  the  results  of  sand  box  experiments  to  "confirm"  his 
results  for  a  cross  section  of  a  typical  Foothills  structure — his  Alpha  structure. 

The  sand  box  is  an  unsatisfactory  piece  of  apparatus.  Haddow  (1965) 
pointed  out  that  the  form  of  the  "fault"  surfaces  developing  in  the  sand,  assuming  the 
sand  behaves  as  a  perfectly  plastic  body  (Ode,  1960)  is  strongly  dependent  on  the  shape 
and  frictional  properties  of  the  surface  of  the  plunger.  In  particular  if  the  plunger 
is  smooth  and  planar,  the  contact  of  the  plunger  with  the  base  of  the  sand  box  is  a 
singular  point  in  a  plastic  slip  line  field.  This  field  has  at  least  one  simple  solution, 
a  fault  plane  at  45°  to  the  horizontal  originating  at  the  singular  point  and  reaching 
the  free  surface  of  the  sand.  This  solution  strongly  resembles  Link's  illustrated  models, 
but  since  there  is  no  equivalent  of  the  plunger  in  the  Rocky  Mountain  orogeny  it 
hardly  confirms  Link's  interpretations. 

A  more  reasonable  approach  was  that  of  Hafner  (1951)  who  attempted  to 
explain  some  of  Hubbert's  (1951)  results  with  sand  box  experiments  using  the  theory 
of  elasticity.  Hafner  (p.383)  visualized  the  following  set  of  events  in  the  failure 
of  a  rock  body  under  a  stress  field. 

A  condition  of  stress  is  produced.  As  long  as  the  shearing  stress  is  every¬ 
where  below  the  elastic  limit  of  the  material,  the  stress  field  can  be  determined. 
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With  further  increase  of  the  stresses,  the  yield  strength  is  exceeded  at  some  point  and 
faulting  occurs.  The  positon  of  this  point  can  be  determined  from  the  initial  condit¬ 
ions.  Around  the  point  there  is  a  drastic  reduction  of  stress,  while  further  away  the 
stress  field  remains  substantially  unchanged.  The  fault  thus  provides  local  relief 
in  the  most  intensely  stressed  portion  of  the  body.  The  further  build-up  of  the  stress 
field  is  presumed  to  go  on  until  renewed  faulting  eventually  takes  place  at  some  distance 
from  the  first  fracture  where  the  stress  relief  was  negligible.  So  the  new  fracture  still 
conforms  quite  closely  to  the  old  stress  field. 

The  opinions  set  forth  in  Hafner's  paper  have  prompted  the  use  of  his 
Fig.  6  to  give  a  qualitative  explanation  of  the  stress  field  in  the  Rockies  at  the  time 
of  the  orogeny  (Charlesworth  1959,  Fig.  4).  Fig.  6b  shows  that  failure  takes  place 
first  in  the  top  left  hand  corner  of  the  block.  If  the  effect  of  the  failure  is  local, 
faulting  resumes  at  other  localities  on  the  same  "thrust  fault  surface"  (Fig.  6b) 
outside  the  range  of  the  local  stress  relief.  There  seems  no  reason  why  this  process 
should  stop  till  the  block  is  actually  split  along  this  trajectory.  The  stress  this  block 
can  sustain  is  then  simply  the  sliding  friction  on  the  fault  plane.  Stresses  greater  than 
this  amount  cause  movements  of  the  block  and  Hubbert  in  the  immediately  preceding 
paper  has  shown  that  commonly  experienced  movements  are  caused  by  negligible 
increases  in  the  stresses.  There  is  then,  no  reason  why  the  stress  in  the  region  of  the 
fault  should  again  approach  the  value  causing  initial  failure  in  the  rock  nor  will 
stresses  capable  of  causing  failure  be  transmitted  into  the  stable  mass  to  the  right  so 
the  situation  which  Hafner  illustrates  in  Fig.  6b  is  improbable.  At  x  =  -4  the  rock 
mass  has  failed,  at  x  =  -12  the  already  fractured  rock  mass  is  supporting  a  stress  three 
times  its  strength.  If  c  (Fig.  6b)  is  reduced  to  a  more  likely  value,  the  boundary 
between  the  stable  and  the  unstable  fields  approaches  the  horizontal  and  the  thrust 
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planes  lose  their  curvature.  So  much  of  the  original  close  resemblance  to  the  struct¬ 
ural  cross  sections  of  the  Rockies  is  lost.  As  Hafner  himself  concludes  . . .  ."the  examples 
are  only  rigorously  correct  if  the  stress  is  everywhere  below  the  strength  of  the  material." 

Attempting  to  use  elasticity  theory  to  analyze  the  behaviour  of  the  rock 
mass  has  several  disadvantages.  At  the  root  of  these  is  the  inability  of  the  theory 
to  deal  with  bodies  in  which  finite  displacement  takes  place  under  stresses  above  the 
elastic  limit. 


Above  the  elastic  limit  of  the  substance  some  aspects  of  the  body's  be¬ 
haviour  might  be  reasonably  expected  to  be  approximated  by  a  perfectly  plastic  body. 
The  model  to  be  examined  will  be  a  perfectly  plastic  body  undergoing  plane  strain. 

It  is  slipping  across  a  perfectly  rough  horizontal  surface  towards  a  rigid  area.  Much 
of  the  analysis  follows  Jaeger  (1962,  pp.  146-149) 

The  initial  object  of  the  analysis  is  to  determine  the  form  of  potential 
surfaces  of  faulting  in  the  moving  mass.  Ode  (1960)  pointed  out  that  in  plasticity 
analysis,  faults  can  be  interpreted  as  surfaces  of  velocity  discontinuity  along  which 
movement  has  taken  place  for  some  time. 

1)  The  principal  axes  of  rate  of  strain  and  of  stress  deviation  are  the  same. 

2)  The  substance  is  uncompressible, 

e  +-  e  +  €  =0 

I  2-  2> 

where  the  are  the  principal  strains. 

3)  Elastic  strains  are  negligible  in  comparison  with  the  plastic  strains. 

4)  The  principal  stress  deviations  are  proportional  to  the  principal  rates  of  strain. 


That  is 


S  =  2cpe 


where  the  S.  are  the  principal  components  of  the  stress  deviation  at  some  fixed 
i 

point,  the  e.;are  the  principal  rates  of  strain  at  the  point,  cj)  is  a  constant. 
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The  above  are  Jaeger  s  (Jaeger,  1962)  ^  29  (3).  He  shows  that  they  lead  to: 

^  =  2  fa  ,  Sy  =  zfa  >  sz  =  2  cj>i2 


Jaeger  S  29  (4). 


^  Vy* >  ^ 


rj 


Y> 


r 


In  plane  perfectly  plastic  strain  the  z  direction  is  conventionally  chosen 
as  the  direction  of  zero  motion  . 


cr  r 

z 


from  Jaeger  S  29  (4) 


S_  -  O 


As  the  stress  deviator  in  this  direction  is  zero,  the  z  direction  must  be 


a  direction  of  principal  stress 


By  Jaeger  S  27  (4) 


<f  =  <y  +  cr' 

3  I _  2. 


CT  >  C T  are  the  principal 

stresses  in  the  xy  plane, 


Besides  the  assumptions  involved  in  the  concept  of  perfect  plasticity 
some  yield  condition  must  be  assumed  before  the  analysis  can  proceed.  This  yield 
condition  holds  at  every  point  wi  thin  the  region  in  which  the  body  is  behaving  per¬ 
fectly  platically.  The  most  common  choice  of  yield  condition  is  that  due  to  Von  Mises:- 

J^  =  a  constant, 

2  2  2 

J  =  (S.  +  S0  +  S*  ),  J«  is  the  second  invariant  of  the  stress 

2  1  2  o  4 

deviator  tensor, 
by  Jaeger  ji  27  (24) 

T 

where  is  the  yield  stress  in  uniaxial  tension 

(or  compression). 

By  Jaeger  $27  (25)  and  $27  (4),  in  plane  perfectly  plastic  strain, 
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where  k  =  CT"  •  3  ^ 

Since  the  maximum  value  of  the  shear  stress  is 

i  ((°*  ~°'y)  ^  ^T7y)  (2) 

with  the  conventions  of  sign  and  definitions  given  by  Jaeger  (1962),  pp.  3-4) 
(1)  becomes  (3) 

There  are  two  other  equations  in 
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*  N  y  ,  •  xy 
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- X?  = 0 


9y 


(4) 


(5) 


X,  Y  are  the  body  forces  per  unit  volume  in  the  x,  y  directions,^)  is  the  density  of 
the  body.  (4)  and  (5)  are  the  equivalents  in  plain  strain  of  Jaeger's  34  (1)  and  (2) 
(Jaeger,  1962).  Jaeger  derives  these  equations.  They  are  equivalent  to  the  assump¬ 
tion  that  the  velocities  U,  V  in  the  x,  y  directions  at  a  point  in  the  plane  do  not  change 
with  time — in  other  words,  the  system  of  stresses  is  in  equilibrium. 

Since  it  is  only  intended  to  consider  gravity  here 
Y  =  O,  X  =  -g 

If  (4),  (5)  are  partially  differentiated  with  respect  to  y,  x  (6),  (7),  result: 
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(8)  results  by  subtraction  and  using  (3) 
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gv~/- 
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<?x  9y 


?\/  (8) 


,nd 
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(8)  is  a  2  order  partial  differential  equation  in  T  .  Normally  its  solution  presents 
considerable  difficulty  but  in  this  case  a  plausible  simplifying  assumption  can  be  made. 
Choosing  x,  y  axes  so  that  the  perfectly  rough  surface  is  given  by  y  -  -a,  the  maximum 
shearing  stress  will  be  developed  along  this  surface~“as  it  has  been  already  assumed 
that  this  surface  is  one  on  which  slip  is  taking  place.  Then  from  the  discussion  of 


Jaeger  S  27  (25), 


/ 


a  boundary  condition  independent  of  x 
of  y  only.  If  this  is  the  case, 

(8)  then  gives  (9) 


This  suggests  that  T  may  be  a  function 


7jL*y  =  0  (9) 

J  =  Ay  +  B,  as  y  =  O  is  the  earth's  surface 
xy 

B  =  O  since  the  atmosphere  cannot  withstand  a  shearing  stress. 

Substituting  the  boundary  conditions, 

T  =  k  when  y  -  -a,  if  the  material  is  slipping  to  the  right  across  the  rough  surface, 
xy 


(10) 
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can  be  found  by  substitution  in  (4)  and  (5) 
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o  y  -  fay  -+  j  (x )  +  c  ('2) 

where  f(y),  g(y)  are  functions  of  x,  y  to  be  determined  and  c  is  a  constant,  Substit- 
uting  (1 1)  and  (12)  in  (3) 


By  comparing  coefficients  of  x,  y. 


$(*)  =  _]£  x  os) 

a 

Substituting  (14)  and  (15)  in  (11)  and  (12) 


In  this  analysis  the  negative  root  is  chosen  in  (18).  This  choice  will  be  justified  when 
the  velocity  field  is  examined. 

Having  derived  the  stresses,  the  next  step  is  to  derive  the  relationship 
between  these  and  the  trajectories  of  maximum  shear  stress.  These  will  coincide  with 
the  directions  of  maximum  shear  strain  and  (according  to  Jaeger  1962,  p,  144)  be 
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the  locus  of  points  at  which  the  relative  movement  between  neighbouring  parallel 
points  is  a  maximum.  They  are  the  "faults"  or  slip-lines  whose  form  was  to  be  derived 
The  trajectories  of  maximum  shearing  stress  bisect  the  principal  axes  of 
stress.  The  equations  connecting  the  stresses  across  an  element  inclined  at  0(~  to 
the  principal  axes  with  the  principal  stresses  are  (Jaeger,  pp.  6-7) 
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The  trajectories  of  maximum  shearing  stress  will  be  inclined  at 
x  axis  where  Ji  =  o(  +  ~^/  (Jaeger,  Fig.  52). 
Replacing  0(  by  j ^  -  ~7T/  ^ 

—  c/  by  k  (from  (1^(22),  (23),  (24)  result. 
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the  gradient  of  one  slip  line  system,  then  from  (22-24) 
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the  other  slip  line  system  is  orthogonal  to  this,  its  gradient 
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So  the  stresses  may  now  be  substituted  into  (^5)  or  (26)  to  derive  the  differential 
equations  of  the  systems  of  slip  lines.  Using  (26) 
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(27) 


dx  (I  -  ) 

Using  the  substitution,  y  =  (X  COS©  the  integral  becomes 


X  =  —  c\ 


(J  -  CoS  6)  d  B 


and  hence  'y  ~  a  cos  X  =  (28) 

(where  D  is  a  constant)  is  the  solution. 

(28)  are  the  parametric  equations  of  a  family  of  cycloids.  From  (27), 

±L  is  always  positive,  so  only  the  portions  where  y  increases  with  x  may  be  taken, 
dx 

The  other  slip  line  system  is  orthogonal  to  these,  it  has  parametric 

equations 

y  =  c\  CoS  b }  x  =  ^  ■+  S“fr\  6)  +■  t  (29) 

E  is  a  constant . 

A  velocity  distribution  compatible  with  the  stress,  (16)  -  (18),  must  now 


be  chosen . 


The  components  of  the  stress  deviator  are  by  (16)  -  (18) 
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These  are  connected  to  the  velocities  at  (x,  y)  by 
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where 


<»  % 

£>  6 

*  ) 


y  are  the  rates  of  extension  in  the  x,  y  directions  and 
is  the  rate  of  shear. 


Y  v 

A  velocity  distribution  compatible  with  compression  in  the  x  direction  i; 
required.  Since  W  ^  I  -  )  is  always  positive,  by  (30),  the  negative 

root  must  be  chosen  in  (16)  -  (18).  From  (31) 


3x 
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(34) 


where  U,  V  are  the  velocities  in  the  x,  y  directions. 

P  is  an  unknown  function  involving  the  strain  rate.  The  choice  of  P  will 
give  some  freedom  to  the  choice  of  velocity  distribution.  There  is  one  boundary 


condition — at  y  =  -a,  V  -  O. 

V  =  xY.  +  N/o 

C\ 


(35) 


where  is  the  upward  velocity  of  the 
free  surface,  satisfies  this  condition  if 


we  choose 

P 

Substituting  (36)  in  (32) 
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(34)  determines  f(y)  , 
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F  is  a  constant. 
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U  =  -2V,(±  +  (l  -%lf)  +1 -  (3') 

The  only  other  analysis  of  this  type  is  the  geological  literature  which  could 
be  found  (Varnes,  1962),  ended  with  the  calculation  of  the  slip  line  fields,  and  did 
not  go  on  to  consider  the  effects  of  finite  displacements  along  slip  lines. 

Fig.  II  illustrates  the  form  of  the  slip  line  field  calculated  in  this  model 
and  represents  diagrammatical ly  the  distribution  of  stresses  on  horizontal  and  vertical 
sections  through  the  body. 

The  "thrust  fault"  in  the  slip  line  field  bears  a  close  resemblance  to  a  typica 
Foothills  thrust — it  rises  very  gradually  out  of  the  zone  of  decollement  and  then  steepens 
more  rapidly  towards  the  surface.  This  superficial  resemblance  is  obviously  a  prelim¬ 
inary  requirement.  At  some  time  in  the  future  the  form  of  a  typical  Foothills  thrust 
should  be  compared  statistically  to  the  form  of  the  ideal  slip  line  and  of  any  other 
models  that  may  have  been  proposed  by  then. 

Hill  (1950)  showed  that  the  calculations  involved  (Chapter  VI-7)  in 
determining  the  velocity  field  around  the  velocity  discontinuity  are  not  so  much  theor¬ 
etically  formidable  as  arthmetical ly  tedious. 
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APPENDIX  1 

THE  ANALYSIS  OF  THE  PRECISION  PARAMETER,  K 

It  has  been  suggested  (on  p.  9  )  that  it  would  be  convenient  to  have  some  means 
of  estimating  the  precision  parameter,  K,  at  a  particular  outcrop  before  visiting  it. 
With  this  in  mind  K  values  for  108  outcrops  where  samples  of  ten  or  more  readings 
had  been  taken  were  analyzed.  (The  restriction  on  the  sample  size  follows  from 
the  statistical  assumption  that  the  error  in  the  determination  of  the  independent 
variable  should  be  negligible).  A  plot  of  part  of  the  data  suggested  the  mean  dip  and 
the  lithology  were  important  in  determining  K  at  a  particular  outcrop. 

To  test  this  suggestion  Dr.  K.  W.  Smillie  proposed  that  for  each  lithology  an 
attempt  be  made  to  fit  the  same  function  to  the  K  values  plotted  against  dip  which 
could  be  transformed  to  a  straight  line,  and  then  examine  the  straight  lines  to  deter¬ 
mine  whether  they  were  coincident.  The  first  step  was  designed  to  show  how  the  value 
of  K  depended  on  the  dip  and  the  second  to  test  whether  the  lithology  had  any  signif¬ 
icant  effect. 

A  preliminary  plot  suggested  the  functions 


x2  =  ax,b  (1) 

X2  =  A  +-®  (2) 

X1 

X2  =  A  +  Blog1QX1  (3) 

might  be  suitable.  These  transform  to 

where 

(1)  x5=  a'  +  b,x4 

(2)  X2  =  a'  +  b'x3 

(3)  X2  =  a'  +b'x4 


)  X2=  K 

) 

)  Xj  =  Dip 

)  A,  B  are  constants 


Xc  =  log.nX„  A  ,  B  are  constants 
o  1U  / 

X4  =  lo910Xl 
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Dr.  Smillie's  multiple  regression  program  (1965,  Chapter  6)  was  used  to  test  the 
the  fit  of  these  functions  to  the  data.  The  results  are  set  out  in  the  Table  1 :2. 

The  results  of  fitting  the  first  function  to  each  of  the  five  lithologies  were  analyzed  by 
the  technique  set  out  in  Hald  (1952),  pp.  579“584.  The  calculations  involved  are 
lengthy  and  only  the  results  will  be  quoted  here.  The  data  for  the  calculations  are 
set  according  to  Hold's  Table  18:19  with  his  notation  in  Table  1 :1  . 

While  the  slopes  of  the  regressions  did  not  significantly  differ,  the  statistics 
showed  at  the  1%  significance  level  that  the  regression  lines  did  not  coincide.  It 
appears  then  that  the  lithology  of  the  strata  has  a  significant  effect  on  the  value  of 
K  as  does  the  dip. 

Suppose  that  for  each  lithology,  the  number  of  samples  of  bedding  planes 
of  this  lithology  is  and  log^  (mean  dip)  of  these  samples  takes  values  x. 

(i=l , . ,  )  and  I og ^ q  K  takes  values  y.  (i=l , . .  . ,  N^)  then.  Hold's  notation  is 


that 
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In  Table  1 :2 
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R  is  the  correlation  coefficient 
F  is  the  critical  F  value 
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TABLE  1:1 


Lithology 

W 

Fernie 

5.365 

12.857 

0.994 

Mount  Head  & 
Turner  Valley 

9.911 

23.644 

1 .558 

Pekisko 

36 . 754 

70.156 

3.914 

Banff 

18.926 

35.888 

0.908 

Pall  iser 

77.105 

148.327 

5.646 

SPb 

V/,  k| 

f 

M 

-2.040 

-2.049 

0.522 

3 

0.1741 

5 

-0.248 

-0.537 

0.263 

7 

0.0376 

9 

-1 .888 

-0.480 

1 .653 

24 

0.0375 

26 

-0.653 

-0.717 

1.057 

11 

0.0961 

13 

-4.406 

-0 . 794 

3.489 

53 

0.0658 

55 

TABLE  1:2 


Lithology 

F 

A  1 

B1 

Fernie 

88.89 

24.00** 

1 

4.771 

-2.049 

47.55 

2.72 

2 

60.69 

4.63* 

3 

Mount  Head  & 

33.54 

3.53* 

1 

3.219 

-0.537 

Turner  Valley 

17.88 

1.23 

2 

25.61 

2.41 

3 

Pekisko 

35.31 

13.10** 

1 

3.377 

-0.480 

31.33 

10.95** 

2 

34.85 

12.84** 

3 

Banff 

30.65 

4.86* 

1 

3.804 

-0.717 

10.82 

1.13 

2 

17.38 

2.31 

3 

Pa  1 1  iser 

50.07 

53.15** 

1 

3.810 

-0 . 794 

29.62 

22.31** 

2 

31.16 

23.98** 

3 

**  significant  at  the  1% 

level 
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TABLE  1:3  DATA  USED  IN  ANALYSIS  OF  K 


EXPOSURE 

DIP  OF 

DIRECTION 

PRECISION 

SAMPLE 

NUMBER 

BEDDING 

OF  DIP 

PARAMETER 

SIZE 

Lithology, 

Palliser-limestone, 

dolomite 

068* 

59.3 

058.6 

0093.5 

16 

071 

10  .6 

234.4 

0524.9 

16 

210 

16.4 

046.9 

3301.4 

11 

254 

44.8 

053.6 

0193.1 

10 

260 

21.5 

176.1 

0670.0 

10 

285 

10.7 

215.8 

0586.6 

10 

291 

46.9 

027.2 

0247.9 

10 

301 

45.9 

040.8 

0498.6 

10 

302 

81.1 

215.2 

0219.0 

10 

303 

75.4 

037.4 

0137.1 

10 

304 

76.8 

229.8 

0189.0 

10 

310 

60.9 

045.2 

0274 . 6 

10 

311 

69.0 

221.6 

0137.6 

10 

312 

49.1 

047.4 

0209.5 

10 

314 

44.5 

221.6 

0122.5 

10 

315 

65.6 

045.4 

0858 . 8 

10 

316* 

68.5 

217.6 

0171.4 

10 

317* 

61.2 

213.6 

0341  A 

10 

318 

59.3 

039.6 

0129.5 

10 

319 

37.7 

216.9 

0141.9 

10 

320 

61.4 

036.3 

0393.4 

10 

321 

44.6 

037.4 

0385.6 

10 

323 

64.9 

026.1 

0047.3 

10 

330 

06.9 

186.6 

0826.4 

10 

331 

10.2 

239.1 

0700 . 1 

10 

339 

14.4 

205.7 

0644 . 3 

10 

340 

13.0 

211.5 

0973.1 

10 

341 

14.3 

215.1 

1840.9 

10 

342 

21.7 

217.9 

1409.1 

10 
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TABLE  1  :3  continued 


EXPOSURE 

NUMBER 

DIP  OF 
BEDDING 

DIRECTION 

OF  DIP 

PRECISION 

PARAMETER 

SAMPLE 

SIZE 

344 

85.1 

221.1 

0532.2 

10 

344 

15.5 

216.0 

0405.1 

10 

345 

50.9 

225.5 

0109.0 

10 

345 

11.8 

206.4 

1857.9 

10 

346 

12.3 

206.1 

0960.4 

10 

346 

16.9 

047.7 

0468 . 1 

10 

347 

43  o  6 

043.3 

0486 . 7 

10 

347 

10.0 

211.9 

0589.8 

10 

349 

11.3 

212.1 

0784.5 

10 

350 

12.6 

218.0 

1408.4 

11 

351 

17.6 

223.2 

0521.9 

10 

352 

18.6 

217.7 

0565.6 

10 

356 

10.0 

233.3 

0755.1 

11 

359 

71.6 

026.1 

0309  7 

10 

360 

65.1 

032.1 

0303.0 

10 

364 

20.5 

034.6 

0234 . 9 

10 

374 

26.3 

019.7 

1082.3 

10 

375 

25.8 

021.8 

0407 . 5 

10 

376 

11.6 

009.7 

1131.8 

10 

378 

18.1 

029.1 

1237.0 

10 

381 

18.8 

005.9 

0480 . 8 

10 

382 

14.9 

026.7 

0825.5 

10 

383 

23.1 

020.8 

1104.1 

10 

388 

09.5 

225.6 

1461.3 

11 

390 

81.7 

027.9 

0137.2 

10 

391 

30.6 

022.6 

0258.9 

10 

401 

23.1 

184.3 

0249.4 

10 

402 

14.6 

190.3 

0659.9 

10 

404 

07.7 

185.7 

2450 . 2 

10 

405 

18.1 

207.2 

0576.0 

10 
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EXPOSURE 

DIP  OF 

DIRECTION 

PRECISION 

SAMPLE 

NUMBER 

BEDDING 

OF  DIP 

PARAMETER 

SIZE 

Lithology, 

Banff-argil  laceous 

limestone,  shale 

066* 

58.2 

214.8 

0229.3 

10 

no 

77.9 

038.4 

0392.6 

11 

111 

78.4 

036.1 

0423.9 

10 

124 

30.6 

055.8 

0562 . 3 

12 

274 

15.4 

214.6 

1203.8 

10 

276 

14.0 

216.9 

0594 . 2 

10 

277 

23.1 

209.1 

0935.8 

10 

278 

24.6 

209.7 

1860.9 

10 

296 

70.4 

021.5 

0111.6 

10 

300 

43.3 

058.3 

0259.0 

10 

332 

21.3 

222.1 

0252.4 

10 

333 

13.9 

212.1 

0606 . 3 

10 

407 

19.1 

203.8 

0726 . 3 

10 

408 

21.6 

205.0 

2060 . 9 

10 

Lithology, 

Pekisko- limestone 

042 

15.2 

027.5 

0236.3 

67 

046 

08.4 

010.2 

0686 . 4 

10 

101 

07.4 

020.6 

1555.3 

15 

104* 

86.5 

214.4 

0159.9 

27 

107 

15.4 

018.2 

1017.7 

16 

115 

18.6 

299.0 

1363.3 

11 

116 

16.1 

212.2 

0804.2 

12 

120 

76.6 

037.5 

0370.2 

15 

121 

74.3 

039.4 

0599.5 

16 

123 

79.5 

043.7 

0561 .9 

20 

262 

22.0 

188.2 

0331.0 

10 

267 

10.2 

231.7 

0926 . 3 

14 
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TABLE  1:3  continued 

EXPOSURE 

DIP  OF 

DIRECTION 

PRECISION 

SAMPLE 

NUMBER 

BEDDING 

OF  DIP 

PARAMETER 

SIZE 

269 

10.4 

234.3 

1258.7 

12 

271 

15.2 

206.3 

1654.9 

10 

273 

32.8 

208.8 

0544.7 

10 

279 

22.3 

200.5 

0657.3 

10 

289 

54.8 

035.7 

0248.9 

10 

299 

71.8 

047.9 

0163.8 

10 

308 

70.3 

048.0 

0209.7 

10 

335 

08.9 

206.7 

0748 . 1 

10 

336 

13.8 

198.7 

0872.4 

10 

337 

10.8 

198.9 

172.7 

10 

338 

08.5 

215.9 

483.4 

10 

363 

78.0 

036.3 

248.1 

10 

365 

83.9 

029.0 

297.2 

10 

387 

73.5 

018.5 

421.1 

10 

307 

81.7 

043.8 

124.6 

10 

Lithology, 

Turner  Valley,  Mount 

Head  -  vuggy, 

massive  bedded  dolomite 

003 

08.6 

038.5 

216.1 

14 

004 

05.0 

035.3 

569.1 

11 

008 

08.2 

030.4 

467.4 

14 

027 

21.1 

006.5 

241.7 

19 

057 

07.1 

051.7 

1033.6 

12 

060 

12.6 

003.5 

613.7 

11 

118 

23.4 

039.3 

189.5 

15 

170* 

67.8 

207.3 

2411.5 

12 

255* 

77.7 

213.9 

504.0 

10 

256* 

77.8 

212.5 

1129.7 

10 

258 

16.4 

208.1 

460.8 

12 

261 

24.8 

188.7 

395.2 

10 
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TABLE1 :3  continued 


EXPOSURE 

DIP  OF 

DIRECTION 

PRECISION 

SAMPLE 

NUMBER 

BEDDING 

OF  DIP 

PARAMETER 

SIZE 

263 

84.8 

222.0 

209.3 

10 

286 

61,8 

204.6 

123.9 

10 

344* 

85.1 

221.1 

523.2 

10 

366* 

66 .6 

208.5 

176.7 

10 

Lithology, 

Fernie-Conglomerate, 

mudstone,  chert, 

sandstone 

049 

04,2 

062.2 

1919,6 

10 

050 

04.6 

313.4 

1247.0 

10 

0  77 

09.2 

242.4 

2364 . 1 

17 

169* 

54.2 

194.8 

680.6 

17 

255* 

74.8 

217.8 

278.6 

10 

256* 

74.8 

217.8 

278.6 

10 

259 

18.8 

216.4 

212.2 

16 

263* 

81.3 

223,3 

2689.7 

10 

288 

68.0 

207,3 

6.3 

10 

298* 

88.5 

227.9 

424.7 

10 

305* 

86,2 

043,5 

462.8 

10 

343* 

83.3 

222.4 

1402.6 

10 
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APPENDIX  2 

TO  SHOW  THAT  THE  PROBABILITY  OF  THE  INTERSECTION  OF  TWO  PLANES  THE 
NORMALS  TO  WHICH  HAVE  FISHER  DISTRIBUTIONS  DOES  NOT  ITSELF  ALWAYS 
HAVE  A  FISHER  DISTRIBUTION. 

It  is  sufficient  to  cite  a  particular  case  in  which  the  normals  do  not  have 
a  Fisher  distribution . 

Consider  two  coincident  planes.  Then  the  intersection  of  the  planes  is 
equally  likely  to  lie  in  any  direction  within  the  two  planes.  If  the  normals  to  the 
planes  have  Fisher  probability  distributions,  this  probability  will  be  distributed  with 
axial  symmetry  around  the  normals.  The  probability  of  the  intersection  of  the  planes 
lies  in  a  band  symmetrical  about  the  planes  and  with  a  maximum  density  in  the 
planes.  This  distribution  is  illustrated  in  Ramsay  (1964)  Figs.  3-5. 

While  the  distribution  of  the  probability  of  the  intersection  thus  has  axial 
symmetry  about  the  normal  to  the  planes  this  coincides  with  a  minimum  likelihood  of 
the  intersection.  The  distribution  of  the  intersection  is,  therefore,  not  a  Fisher 


distribution . 
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APPENDIX  3 

THE  MISSISSIPPI  AN- JURASSIC  UNCONFORMITY  NEAR  NORDEGG,  ALBERTA 

INTRODUCTION 

In  the  Rocky  Mountain  Foothills  of  central  Alberta,  the  essentially 
concordant  contact  between  the  Mississippian  Mount  Head  and  Jurassic  Fernie  Form¬ 
ations  is  exposed  near  Nordegg  in  a  low  cut  on  the  south  side  of  Highway  1 1 ,  quarter 
of  a  mile  east  of  the  bridge  over  the  abandoned  railway  line  in  Shunda  Gap.  This 
exposure,  in  which  the  beds  are  slightly  overturned,  was  examined  to  determine  the 
angular  relationship  between  the  two  formations. 

DETERMINING  ANGULAR  RELATIONS  BETWEEN  SUB-PARALLEL  SURFACES 

Determining  the  angular  relations  between  two  planar  structures  is  diff¬ 
icult  if,  as  in  the  case  of  bedding  in  the  Mount  Head  and  Fernie  Formations  at  the 
above  road  cut,  the  structures  are  almost  parallel.  Under  these  conditions  it  is  diff¬ 
icult  to  decide  whether  an  observed  difference  in  orientation  is  significant,  or  whether 
it  results  from  scatter  caused  by  a  combination  of,  say,  instrument  error,  observer 
error,  and  the  departure  of  the  bedding  planes  from  perfect  planarity.  The  problem 
of  how  to  recognise  significant  differences  can  be  approached  by  putting  forward  a 
plausible  hypothesis  about  the  distribution  of  this  scatter.  The  hypothesis  is  that 
repeated  measurements  of  a  planar  structure  within  a  domain  (a  rock  body  that  is 
homogeneous  with  respect  to  the  planar  structure  being  studied)  have  a  distribution 
about  the  estimated  mean  of  the  measurements  which  is  given  by  the  Fisher  (1953) 

probabi I ity  density  function ,  P  (©■)/ 

P  (9)  =  BeKc0s6, 
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where 

B  is  a  constant, 

9  is  the  angle  between  a  given  measurement  and  the  estimated  mean,  and 

K  is  the  precision  parameter,  the  precision  parameter  is  the  three-dimen¬ 
sional  analogue  of  the  inverse  of  the  variance,  such  that  the  lower  the  dispersion  of  a 
sample  of  measurements  the  higher  its  precision  parameter. 

Three  domains  elsewhere  in  the  Nordegg  area  were  investigated  with  a 
view  to  testing  the  validity  of  the  hypothesis.  In  each  the  distribution  of  the  60  to  70 
measurements  of  bedding  does  not  differ  significantly  from  a  Fisher  distribution  .  It 
is  therefore  assumed  that  any  number  of  randomly  collected  measurements  of  a  planar 
structure  in  any  domain  have  a  Fisher  distribution. 

If  the  two  samples  to  be  compared  have  Fisher  distributions  it  is  poss¬ 
ible  to  set  up  an  approximate  F-test  based  on  an  analysis  of  variance  analogy.  This 
is  the  test  of  the  null  hypothesis  that  the  two  samples  come  from  a  population  with  the 
same  mean,  i.e.  that  their  orientations  do  not  differ  significantly.  The  hypothesis 
is  conventionally  rejected  if  there  is  less  than  once  chance  in  twenty  of  its  being 
true.  The  rejection  of  the  null  hypothesis  implies  that  the  two  samples  come  from 
populations  with  different  means  and  hence  that  their  orientations  are  significantly 
different.  According  to  Watson  (1956),  for  samples  with  Fisher  distributions  and 
equal  K  values,  the  hypothesis  is  rejected  if 
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(l^i7  nij.  ,  n  ^ .)  and  (^.z  it^.z  n2 ^>e' n9  the  direction  cosines  of  the  normals  to 
the  planar  structure  (upward  directed  in  the  case  of  bedding)  to  the  i^  measurement 
in  samples  1  and  2  respectively, 
and  N  =  +  l\Lz 

N,  and  being  the  numbers  of  measurements  in  samples  1  and  2  respectively, 

2.  /  J±i  Mj.  \2  /  N/,  JLu  \2  ,  JV,  tfjr  v2- 

R  -  (2 i„^iu )*(2V +Z*d 

and  F2  2((sj-2)  's  tabulated  'n  Lindley  and  Miller  (1953,  p.8). 

The  estimated  precision  parameters  K^,  and  K^,  f  k*  — 

V  '  m  — 


N,  -  R,  , 

of  any  two  samples  to  be  compared  are  unlikely  to  be  equal,  which  suggests  that  the 
F-test  described  above  will  usually  be  inapplicable.  However  these  estimates  are 
subject  to  large  errors  (Watson  1956,  p.  154)  and  Watson  and  Irving  (1957,  p.  292) 
have  suggested  an  approximate  test  of  the  hypothesis  that  the  estimated  precision 

parameters  of  the  two  samples  are  estimates  of  the  same  precision  parameter  K^,  (  ^K^). 

^  a.  / 

They  commented  that  "values  of  K^/l^  far  from  unity  strongly  suggest  ^  K^." 

,  l<2  are  the  precision  parameters  of  the  populations  from  which  the  samples  were 

taken.  The  effect  of  departures  from  the  hypothesis  on  the  validity  of  the  F-test  has 

not  been  examined. 


A  Fortran  IV  computer  program  to  calculate  the  value  of  the  statistic 
^  )  C  ^  1  ^2. _ ^  (here  called  the  coplanarity  statistic) 

(M  -  f?,  -^) 

from  an  input  of  strikes  and  dips  of  the  planar  structure  is  available.  If  this  value 
is  greater  than  the  corresponding  F  value,  the  null  hypothesis  may  be  rejected  at  the 
conventional  level  of  significance. 

The  angular  discordance  0  between  the  two  samples  can  be  determined  from 
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where  (1^ ,  m^ ,  n^)  and  (l^,  m^,  n^)  are  the  direction  cosines  of  the  estimated 
means  of  the  two  samples  of  the  planar  structure  (upward  directed  in  the  case  of 
bedding),  given  in  the  case  of  sample  1  by 

1  m  .  -  Jfi'N  n,  =  2  R, 


U 


R 


R  .  R 

A  Fortran  IV  program  has  been  designed  to  calculate  estimated  means  of  samples 
(Muecke,  1965). 


THE  UNCONFORMITY  IN  THE  SHUNDA  GAP  ROAD  CUT 


It  can  be  verified  by  substitution  in  the  test  that  for  domains  with  precision 
parameters  of  500,  10  measurements  are  required  in  each  sample  to  distinguish  est¬ 
imated  means  two  degrees  apart.  Previous  experience  had  shown  that  the  precision 
parameters  of  the  samples  would  likely  be  between  200  and  700,  so  10  measurements 
of  bedding  were  taken  on  either  side  of  the  contact  in  the  Shunda  Gap  road  cut. 

The  bedding  planes  on  which  these  measurements  were  made  are  those  closest  to 
(within  15  feet)  the  contact. 

Before  the  test  can  be  applied  it  is  necessary  to  be  confident  that  each 
15-foot-long  outcrop  constitutes  only  one  or  part  of  one  domain,  i.e.  that  each 
sample  of  10  measurements  can  be  considered  to  be  a  random  sample  from  a  single 
Fisher  distribution.  The  precision  parameters  of  two  adjacent  Mount  Head  and 
Fernie  domains  for  which  Fisher  distributions  have  been  established  are  247  and  207, 
respectively  •  Because  both  samples  from  the  road  cut  have  precision 
parameters  of  over  200  (Table  3:1),  it  is  unlikely  that  either  sample  comes  from 
more  than  one  domain.  This  reasoning  is  far  from  rigorous,  but  over  30  readings  are 
needed  in  a  sample  before  its  fit  to  a  Fisher  distribution  can  be  tested  statistically 
(Watson  and  Irving,  1957,  p.  293).  It  was  not  possible  to  collect  samples  this  large 
from  close  to  the  unconformity,  so  the  above  rough  test  had  to  be  applied.  Computer 
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programs  for  the  rigorous  test  are  available. 

The  results  of  the  test  on  the  two  samples  from  the  Shunda  Gap  road  cut 
(locality  1)  are  given  in  Table  3:1  .  It  can  be  seen  from  the  table  that  because  the 
value  for  the  coplanarity  statistic  is  greater  than  the  corresponding  F  value,  the 
hypothesis  that  the  samples  are  from  the  same  population  should  be  refected  at  the 
conventional  level  of  significance,  i.e.  the  Mount  Head  and  Fernie  Formations  at 
Shunda  Gap  are  discordant.  It  should  be  pointed  out  however,  that  because  the 
estimated  precision  parameters  of  the  Mount  Head  and  Fernie  samples  are  markedly 
different,  this  test  may^be  only  approximate. 

The  estimated  mean  orientations  of  bedding  in  the  two  formations  are  given 
in  Table  3:2.  The  probable  orientation  of  the  Mount  Head  at  the  time  of  deposition 
of  the  Fernie,  obtained  by  rotating  the  Mount  Head  bedding  about  an  axis  parallel 
to  the  Fernie  strike  by  an  angle  equal  to  the  dip  of  the  Fernie,  is  in  the  last  column 
in  Table  3:2.  Ramsay  (1  961)  and  Cummins  (1964)  have  discussed  the  conditions  under 
which  such  rotations  are  permissible. 

The  discordant  contact  between  the  two  formations  at  the  road  cut  in 
Shunda  Gap  could  conceivably  be  a  fault  rather  than  an  angular  unconformity.  On 
Shunda  Mountain,  six  miles  northwest  of  this  locality,  a  chert  and  dolomite  pebble- 

* 

conglomerate  rests  on  an  irregular  Mount  Head  surface  (Brady,  1958).  In  the  road 
cut,  two  feet  of  friable,  calcareous  pebbly  sandstone  and  five  feet  of  soft,  thinly 
bedded  mudstone  overlie  thinly  bedded,  buff  dolomites  of  the  Mount  Head  Formation. 
The  Jurassic  succession  then  becomes  increasingly  calcareous  with  the  appearance 
of  thin,  black  limestones.  Therefore  the  contact  between  the  two  formations  in  the 
road  cut  gives  over  appearance  of  being  an  unconformity  and  not  a  fault. 
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THE  UNCONFORMITY  ELSEWHERE  IN  THE  NORDEGG  AREA 

Three  other  localities  where  the  actual  contact  between  Mount  Head  and 
Fernie  beds  is  not  exposed,  but  where  10  measurements  of  bedding  could  be  taken 
on  either  side  of  the  contact  and  within  20  feet  of  it,  were  also  examined.  Local- 
ity  2  is  on  the  old  road  through  Shunda  Gap,  50  yards  west  of  the  railway  bridge 
over  the  road.  Locality  3  is  under  the  railway  bridge  over  Shunda  Creek  on  the 
south  side  of  the  creek.  Locality  4  is  on  the  creek  flowing  out  of  the  gap  between 
Shunda  and  Coliseum  Mountains,  150  yards  south  of  where  the  trail  up  to  the  Shunda 
Mountain  Lookout  crosses  the  creek  for  the  first  time.  The  data  for  localities  2-4 
is  greater  than  the  corresponding  F  value,  here  too  the  attitude  of  the  Fernie  is 
significantly  different  from  that  of  the  underlying  Mount  Head.  The  mean  orient¬ 
ations  of  bedding,  and  the  probable  attitude  of  the  Mount  Head  at  the  time  of  de- 
position  of  the  Fernie,  are  shown  for  each  locality  in  Table  2.  The  similarity  of 
these  data  to  those  for  locality  1  suggest  that  the  contacts  at  localities  2-4  are  also 
unconformities.  The  ratio  of  the  estimated  precision  parameters  at  each  locality 
is  not  far  from  unity,  so  the  F-test  is  probably  applicable  in  each  case. 

CONCLUSIONS 

At  four  localities  in  the  Nordegg  area,  Fernie  strata  appear  to  lie 
with  angular  discordance  on  Mount  Head  beds.  The  angle  of  discordance  ranges  from 
3.2  to  4.2  degrees.  At  the  time  of  deposition  of  the  basal  Fernie  the  average 
orientation  of  the  Mount  Head  as  determined  by  the  vector  mean  of  the  four  det¬ 
erminations  was  approximately  strike  -  N  50°W,  dip  -3.8°SW.  This  conclusion 
agrees  with  that  which  can  be  drawn  from  considerations  regarding  the  overall 
stratigraphy  of  the  area  (see  e.g.  Bokman,  1963,  Fig.  5). 
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In  this  appendix  the  author  puts  forward  a  method  of  recognizing  significant 
small  angular  discordances  between  surfaces.  Bedding-plane  orientation  data  from 
either  side  of  the  Mississippian-Jurassic  unconformity  near  Nordegg,  used  to  i  1 1  us— 
trate  the  method,  point  to  the  angular  nature  of  this  contact.  The  method  might 
also  be  used  to  identify  unconformities. 


ADDITIONAL  REFERENCES 
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Ramsay,  J.G.,  1961,  The  Effects  of  Folding  on  the  Orientation  of  Sedimentation 
Structures,  Jour.  Geol.,  Vol.  69,  pp.  84-99. 
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TABLE  3:1  -  SAMPLES  SIZES  (N^  N^,  R]  and  VALUES  AND  PRECISION  PARAMETER! 

(Kr  K2)  for  mount  head  and  fernie  samples,  r  r  r,  the  coplanarity 

STATISTIC  (C.S.),  AND  THE  CORRESPONDING  F  VALUE. 


Mount  Head  Fernie 


Local  ity 

N, 

Ri 

Ki 

N2 

R2 

K2 

R 

CS  F2,2(N-2) 

1 

10 

9.9570 

209 

10 

9.9967 

2690 

1 9 . 9432 

4.08  3.26 

2 

10 

9.9920 

1130 

10 

9.9933 

1344 

19.9764 

10.90  3.26 

3 

10 

9.9821 

504 

10 

9.9677 

279 

19,9323 

6.27  3.26 

4 

12 

11.9761 

460 

16 

15.9293 

212 

27.8930 

3.50  3.16 
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TABLE  3:  2  -  VALUES  FOR  THE  ESTIMATED  MEAN  ORIENTATIONS  OF  BEDDING 
IN  THE  MOUNT  HEAD  AND  FERNIE  FORMATIONS,  AND  FOR  THE  PROBABLE 
ORIENTATIONS  OF  BEDDING  IN  THE  MOUNT  HEAD  AT  THE  TIME  OF  DEPOS¬ 
ITION  OF  THE  FERNIE,  AT  FOUR  LOCALITIES  IN  THE  NORDEGG  AREA.  OVER¬ 
TURNED  DIPS  ARE  MARKED  WITH  AN  ASTERISK. 


Locality 


1 

2 

3 


Mount  Head  Fernie  Mount  Head  in 

Jurassic  time 


strike 

dip 

strike 

dip 

strike 

dip 

132.0° 

84.8°SW* 

133.3° 

81 ,3°SW* 

85° 

4. 1°$W 

122.5° 

77.8°SW* 

121 .2° 

74. 6°  SW* 

134° 

3.2°SW 

123.9° 

77.7°SW* 

127.8° 

74.8°SW* 

o 

o 

CNI 

4 . 2°SW 

126.4° 

18.8°SW 

118.1° 

16.4°SW 

O 

CO 

3.8°SW 

4 
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APPENDIX  4 


COMPUTER  PROGRAMS 

1)  To  find  the  axes  of  cylindrical  folds  and  to  test  for  their  cylindrical  nature. 

Purpose 

The  purpose  of  these  programs  is  to  delineate  portions  of  a  fold  which  are  cylindrical, 
to  calculate  the  fold  axes  of  such  portions  and  to  place  cones  of  confidence  around  these 
axes . 

Method 

The  statistical  basis  of  the  method  is  derived  in  Watson  (1960)  and  Creer  (1962). 

tH 

Suppose  N.  determinations  of  the  bedding  are  made  in  the  i  of  p  domains  and  yield  a 
mean  normal  with  direction  cosines  (I.,  m.,  n.)  and  length  R..  The  null  hypothesis 
is  that  all  the  p  mean  vectors  lie  in  a  plane.  For  domains  with  equal  precision  parameters 


a  test  of  coplanarity  is  given  by: 


Ni , 


(min .  eigen  value  of  V) 


V  = 


N;-Rj) 

][  R;lfw>|-X*irnia'  W 


The  maximum  likelihood  estimator  of  the  direction  of  the  normal  to  the  best 
fitting  plane  is  the  eigen  vector  corresponding  to  the  least  eigen  value  of  the  Matrix  V. 

The  direction  of  the  axes  of  the  ellipse  of  confidence  round  the  eigen  vector  are 
given  by  the  other  eigen  vectors,  the  minor  axis  corresponding  to  the  eigen  vector  assoc¬ 
iated  with  the  smaller  eigen  value.  The  lengths  of  the  semi-axes  are: 
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where  A  are  the  eigen  values  of  V. 

For  domains  with  unequal  K  (the  more  general  case),  a  test  of  coplanarity  is  given  by 
min.  eigen  value  of  W  < 


V  = 


p  X 

1.  &  ir  2 

2  Bi  n,-  2lBh''Wl  ^ 


J3 


(N;  -e.) 


The  directions  of  the  normal  to  the  best  fitting  plane  is  given  by  the  eigen  vector  assoc¬ 
iated  with  the  least  eigen  value.  The  axes  of  the  ellipse  of  confidence  around  the  eigen 
vector  are  given  by  the  other  eigen  vectors,  the  minor  axis  corresponding  to  the  eigen 
vector  associated  with  the  intermediate  eigen  value.  The  semi  axes  are  given  by; 


A  are  the  eigen  values  of  W. 
Restrictions 


The  input  data  for  this  program  must  be  in  the  following  form: 
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Card  1 

The  number  of  matrices  to  be  processed  is  punched  in  columns  2,3,4. 

Card  2 

t  h 

A  three  figure  number  identifying  the  i  domain  to  the  user  is  punched  in 
2,3,4.  The  number  lies  within  the  range  0-998. 

Card  3  - 

The  program  as  written  considers  only  the  downward  directed  normals  to 
bedding  planes.  When  the  dips  within  a  domain  vary  from  nearly  vertical  to  overtuned 
this  could  lead  to  variations  of  180°  in  the  trend  of  the  normals  and  erroneous  values  of 
R..  Consequently  a  parameter  card  consisting  of  dextral  (punched  in  2,  3,  4)  sinistral 
(in  6,  7,  8)  azimuthal  limits  for  each  i  must  be  read  in.  It  is  sufficient  to  assign  the 
dextral  limit  some  value  between  180  and  360,  and  the  sinistral  limit  some  value  between 
0  and  180,  making  about  90°  with  the  dextral  limit,  then  to  keep  these  values  for  all  i. 
342,  062  are  examples  of  such  limits. 

Data  cards 

The  strike  of  the  bedding  is  punched  into  2,  3,  4. 

The  dip  of  the  bedding  is  punched  into  6,  7. 

A  code  is  punched  into  9. 

The  code  has  the  value  0  when  the  dip  is  90. 

1  when  the  dip  is  to  270-360,  0-89. 

2  when  the  dip  is  to  90-269. 

The  value  of  the  strike  lies  between  0-180. 

The  value  of  the  dip  lies  between  0-90.  No  provision  is  made  for  overturned  rocks. 

tH 

A  strike  of  999  following  a  data  card  signifies  the  card  was  the  N  of  the  N. 
readings  from  a  domain.  It  is  followed  by  card  2  identifying  the  next  station.  This 

is  followed  by  card  3. 
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0 

1 


SIBFTC  MAIN  NOLIST, OECK 

DOUBLE  PRECISION  V(20,20),S(20,20),D(3),PI,RZ,AL,AM,AN,TR,CF,PL, 
- 1 - A IM  »  OL  IMrSLI  M ,  A  ,  6  »  C  ,  E,R,Y,  Z,T  R  6,APLU,PLU  , _ 


4 

21 

FORMAT! IX, 13, IX ) 

3 

NMTRX=0 

6 

111 

C022 J=1 , 3 

7 

00221=1,3 

10 

22 

V( I , J ) =0C0 

13 

CHISQ=3.99COOCOCCOOGCO 

14 

P 1=3. 14159265358979 

15 

RZ=1. CO/180. DO 

16 

CF=PI-*RZ 

17 

A  T  N  =  0  D  0  To  COwp'lje. 

CXY\(k 

20 

21 

istc=o  SoUe  the  matrix  }V 

22 

M  M  T  R  X  =  N  M  T  R  X  + 1 

23 

WRITE(6,217)MMTRX 

24 

217 

FORMAT ( // 1 1H  MATRIX  NO. ,13) 

25 

3  2 

N  =  0 

26 

REAC( 5,41 )NSTN 

30 

41 

FORMAT (  IX  t 13) 

31 

A  L  =  0  D  0 

32 

AM=  0C0 

33 

AN=0C0 

34 

IF  (  NSTIM.LT. 0  )G0  TO  13 

37 

ISTC=ISTC+1 

40 

READ  (5  ,51  )  LL I M , NL I M 

43 

51 

FORMAT ( IX, I 3, IX , 13) 

44 

33 

N  =  N  + 1 

45 

R  E A  C  (5  ,52  )  NSTR  *  NO  I P , NDI R 

51 

52 

FORMAT  (  IX, 13,  IX, 12, IX, ID 

52 

IF(NSTR.EQ.999)G0  TO  11 

, 

55 

PL=90-ND I P 

56 

IF ( NDI R- 1)3, 2, 3 

57 

2 

TR=NSTR+90 

60 

GO  TO  5 

61 

3 

TR=NSTR— 90 

62 

"  I  F  (  T  R  )  4 , 5 , 5 

63 

4 

TR=TR+360 .DO 

64 

5 

CL IM=LL  IM 

65 

IF(TR-DLIM)6,10,10 

6  6 

6 

SLIM=NLIM 

67 

I F ( SL IM—DL IM ) 8, 9,9 

70 

8 

I F ( TR— 90. )10, 10,9 

71 

9 

TR  =  TR+  18C . DO 

72 

PL=PL*CF 

73 

TR=TR*CF 

74 

E=-CS I N ( TR ) 

7  £ 
i 

GO  TO  101  .  _ 

76“ 

'To"" 

"  P  L  =  P  L  *  C  F 

n 

:Q 

TR=TR*CF 

E=DS  IN ( TR ) 

5l 

101 

A=DCOS ( PL ) 

■ 


‘ 
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02 

B  =  D  S  I N ( PL  ) 

03 

C=CCCS ( TR  ) 

04 
\  - 

AL=AL+A«C 

• 

i  5 

AM=AM+E*A  ~ ~ - 

--  -  -  —  ■  ■ 

06 

AN= AN+B 

07 

GO  TO  33 

10 

11 

N  =  N—  1 

11 

TN=N 

12 

R=DSQRT(AL*AL+AM*AM+AN*AN) 

13 

AN=AN/R 

14 

AL=AL/R 

1 5 

A  M = A  M  /  R 

16 

ATN= ATN+TN 

17 

AR=AR+TN-R 

10 

V  (  1  »  2 ) =  V { 1 , 2 ) +R  *AL* AM 

11 

V ( 1 , 1 ) = V ( 1, 1) +R*AL*AL 

12 

V ( 2 , 3 )=V ( 2  , 3 ) +R  * AN* AM 

13 

V( 1,3) = V ( 1,3)+R*AL*AN 

14 

V(2,2)=V(2,2)+R*AM*AM 

15 

V(3,3)=V(3,3)+R*AN*AN 

16 

WRITE  (6, 12 )NSTN 

17 

12 

FORMAT ( / 16H  STATION  NUMBER , 3X , I 3 ) 

50 

WRITE! 6, 15)N 

51 

15 

FORMAT { 17H  NO.  OF  READ INGS , 2X , I 3 ) 

52 

GO  TO  32 

53 

13 

V { 2  ,  1)  =  V( 1,2) 

54 

V ( 3 , 1 ) =V (1,3) 

55 

V ( 3 , 2 ) =V (2,3) 

56 

ARTN=ATN/AR 

57 

WRITE! 6, 14) I STC 

r0 

14 

FORMAT { / 1 7H  NO.  OF  STATI ONS , 2X , I 3 ) 

-1 

CALL  JACOBI (N,V ,S,T0LERC, NORM) 

2 

D03 1CL= 1,2 

3 

C03 10K= 1 , 2 

•4 

IF(DABS(V(K,K) ) . LE. DABS ( V ! K+l , K+l ) ) )GO  TO  310 

7 

003101=1,3 

0 

Z  =  S (K,  I  ) 

1 

S ( K  ,  I  )  =  S ( K+ 1 , I) 

2 

S  (  K  +  1 ,  I  )  =  Z 

3 

Y=V ( K , K ) 

4 

V ( K  ,  K  )  =V ( K+ 1 »  K+l ) 

5 

V { K+ 1 , K+ 1 ) = Y 

. . ......... —  .  — . — - 

6 

310 

CONTINUE 

2 

0066 J= 1 , 2 

3 

00311M=1,3 

4 

nr 

C  (  M  )  =  S  {  J  ,  M  ) 

6 

I F ( C {  1  )  )  36,40,36 

7 

36 

ATRE=D ( 2 ) /D ( 1 ) 

0 

TRE=DATAN ( ATRE) /CF 

1 

GO  TO  37 

2 

40 

I  F  (  C  (  2  )  )  43, 4  2,4  2  ... 

3 

43 

'  T  R  E  =  90 -DO 

4 

GO  TO  420 

5 

4 

—  270  DO  . 

6 

420 

APLU=C ( 3 ) /DSQRT (l.D0-D(3)**2) 

. 
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.77 

PLU  =  C  AT  AN ( APLU) 

o 

o 

PLU=DABS ( PLU/CF ) 

GO  TO  65 

J2 

37 

APLU=D ( 3 ) /DSQRT ( 1 .DO-D ( 3) **2) 

!03 

PLU=DATAN ( APLU) 

'04 

PLU=DABS ( PLU/CF ) 

!  0  5 

I F ( C ( 3  )  )61f60»6C 

'06 

60 

C ( 1  )  =~D ( 1  ) 

07 

C ( 2  )  =-D ( 2 ) 

10 

61 

IF ( D ( 2 ) ) 6 3,62, 62 

11 

62 

I  F ( D  (  1)  )64,65,65 

12 

63 

TRE=TRE+180.D0 

13 

GO  TO  65 

14 

64 

TRE=TRE+360.00 

15 

65 

WRITE(6,214)J,PLU 

16 

214 

FORMAT { / 1 2 H  PLUNGE  OF  ,  I 1 , 8H  EVECTOR , 2X , F7 . 3 ) 

17 

66 

WRITE(6,215)J,TRE 

21 

215 

FORMAT ( 1 1H  TREND  OF  ,I1,8H  E V ECTOR , 2X , F 8 . 3 ) 

22 

STA=ARTN*V (1,1) 

23 

WRITE (6,210) STA 

24 

210 

FORMAT ( /11H  STATIST  IC ,  10X,F10. 4) 

25 

T1  =  DSQRT (CHI  SO/ (  ( V ( 3 , 3 ) -V ( 1 , 1 )  ) » ARTN ) ) 

26 

T2=DSQRT (CHISQ/ ( ( V ( 2 , 2 ) - V ( 1 , 1 ) ) * ARTN ) ) 

27 

T 1 T=T 1 /DSQRT ( 1. D0-Ti**2) 

30 

THET 1  =  DAT  AN (TIT ) /CF 

31 

T2T=T2/DSQRT { 1. D0-T2**2) 

32 

THET2=DATAN(T2T  )/CF 

33 

WRITE(6,212)THET1 

34 

212 

FORMAT ( 2 2 H  RADIUS  OF  CONF I DENCE , 1 X , F 7 • 3 ) 

35 

WRITE(6,212)THET2 

3o 

50 

CONTINUE 

3? 

NMTRX=NMTRX+1 

^0 

IF( NMTRX.LT. MIX) GO  TO  111 

+3  CALL  EXIT 

:»4  END 


■ 


M 

■ 
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0  SIBFTC  MAIN  NOL I  ST , NODECK 

I  DOUBLE  PRECISION  V { 20 , 20 ) , S ( 20 ,20 ) , D ( 3 ) t P I , RZ , AL , AM , AN , TR , CFt PL , 

A TR E, D  L I  M  ,  S  L  I M  ,  A  ,  B  ,  C ,  E , R ,Y, Z, TRE, APLU, PLU , _ 


1 


4 

5 

6 

_ 7 

10 

13 

14 

15 

16 
17 
20 
21 
22 

23 

24 
_26 
27 

30 

31 

32 

35 

36 

41 

42 

43 
47 
50 

53 

54 

55 

56 

57 
60 
61 
62 

63 

64 

65 

66 
67 

70 

71 

72 

73 

74 
J75 

76 

77 
100 
101 


2  CHISQ,TN, T1,T2, THET1, THE T2 , AX, TIT , T2T 

READ(5,21)MIX 

21  FORMAT ( IX  , 13, IX) 

NMTRX=0 

111  D022 J=  1,3 

_ 00221=1,3  _ 

22  V  ( I  »  J  )  =  0  D  0 

1ST  C  =  0  ~ 

MMTRX=NMTRX+1 
WRITE(6,217)MMTRX 

217  FORMAT ( //I 1H  MATRIX  NO. ,13) 

32  N  =  0 


41 


To  Connpi  le,  and 

Solve  iPoL  yyiatrot 

w 


CH  I  SQ= 5. 9900000  0000000 
P 1=3. 14159265358979 
RZ=1. D0/180. DO 
CF=PI*RZ 
READ { 5 ,4  1  )  NSTN 
_Fi 0 R M A T  (  IX, 13) _ 


AL=0D0 
AM=  ODO 
AN=0D0 

IFINSTN.LT. 0)G0  TO  13 
ISTC=ISTC+1 

READ  (5  ,51  )  LLIM,NLIM 


51  FORMAT ( IX, I  3, IX , 13) 

33  N  =  N+ 1 

READ  (5  ,52  )  NSTR , NDI P , NDI R 

52  FORMAT  {  IX, 13, IX, 12, IX, ID 
IF ( NSTR. EQ. 999) GO  TO  11 

_ PL  =  90-NDIP _ _ 

I  F  (  N  D  I  R-  1 )  3 »  2  » 3 

2  TR=NSTR+90 
GO  TO  5 

3  TR=NSTR-90 

I F ( TR ) 4 , 5 , 5 

4  TR=TR+360. DO _ 


5  DL I M=LL  I M 
IF(TR-DLIM)6, 10,10 

6  SL I M  =  NL  I  M 

I F  < SLIM-DLIM)8,9,9 

8  IF( TR-90.D0  )  10,  10,9 

9  _ TR  =  TR+  1 80  .  DO _ 

PL=PL*CF 
TR=TR*CF 
E=— DS I N ( TR ) 

GO  TO  101 

10  PL=PL*CF 

_ TR  =  TR»CF _ 

E=DS  IN ( TR  ) 

101  A=DCOS ( PL ) 

B=DSIN( PL) 

C=DCOS ( TR ) 
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102 

al=al+a*c 

103 

AM=AM  +  E  * A 

~  104 

AN= AN+6 

105 

GO  TO  33 

106 

11 

N=N—  1 

107 

TN=N 

110 

R=DSQRT ( AL*AL+AM*AM+AN*AN ) 

111 

A  N = A  N  /  R 

112 

al=al/r 

113 

AM=AM/R 

114 

AX=R* ( TN-l.DO)/ (TN-R) 

115 

V(3,3)=V(3,3)+AX*AN*AN 

116 

V( 2, 3)=V(2,3)+AX*AN*AM 

117 

V(2,2)=V(2,2)+AX*AM*AM 

120 

V( 1,3)=V( 1 ,3)+AX*AL*AN 

121 

V ( 1 »  2 )=V ( 1 , 2 ) +AX*AM* AL 

122 

V(  1, 1 ) = V  C 1,  1)+AX*AL*AL 

123 

WRITE (6, 12) NSTN 

124 

12 

FORMAT ( / 16H  STATION  NUMBER , 3X • I 3 ) 

125 

WRITEt  6,  15)N 

126 

GO  TO  32 

127 

15 

FORMAT ( 17H  NO.  OF  READ INGS » 2X , I 3 ) 

130 

13 

V ( 2  ? 1 ) = V ( 1,2) 

131 

V ( 3  7 1 ) =V ( 1,3) 

132 

V(3,2)=V(2,3) 

133 

WRITE(6,14) I STC 

134 

14 

FORMAT ( / 1 7H  NO.  OF  ST AT  I ONS , 2X , I 3 ) 

135 

CALL  JACOBI (N,V, S,TOLERC, NORM) 

136 

D03 10L=  1 , 2 

137 

D03 1 0K=  1,2 

140 

IF(DABS(V(K,K)) .LE. DABS ( V ( K+l , K+l ) ) )G0  TO 

310 

143 

D03 101=1,3 

144 

Z  =  S ( K  7 I ) 

145 

S(K,I)=S(K+1,I) 

. 

146 

S ( K  + 1 , I )-Z 

147 

Y=V(K,K) 

150 

V(K»K)=V(K+1,K+1) 

151 

V ( K+l ,K+1  )=Y 

152 

310 

CONTINUE 

156 

157 
160 
162 

163 

164 


311 

36 


0066 J= 1 » 2 
DO 3 1 1M= 1 » 3 
D  (  M  )  =  S  (  J  ,  M ) 

I F ( D ( 1  )  )36,40,36 
ATRE=D( 2)/D( 1 ) 


1  165 

GO  TO  37 

166 

40 

I F ( D ( 2  )  ) 4 3, 42, 42 

167 

43 

TRE=90.D0 

170 

GO  TO  420 

171 

42 

TR£=270.D0 

172 

420 

APLU  =  D ( 3 ) /DSQRT (  1 

173 

174 

175 

176 


37 


PLU=OATAN( APLU) 

PLU  =  DABS( PLU/CF  ) 

GO  TO  65 

APLU=D(3)/DSQRT (l.D0-D(3)**2) 


‘ 


■ 

' 
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177 

PLU=DATAN(APLU) 

200 

PLU=DABS( PLU/CF ) 

201 

I F ( 0(3) )61, 60,60 

202 

60 

D (  1 )  =-D  (  1 ) 

203 

0  (  2  )  =— D  (  2  ) 

204 

61 

IF(D(2))63,62t62 

205 

62 

I F ( 0 (  1 )  )64, 65,65 

206 

63 

TRE=TRE+180.D0 

207 

GO  TO  65 

210 

64 

TRE=TRE+360.D0 

211 

65 

WRITE(6,214)J,PLU 

212 

214 

FORMAT ( /  1 2H  PLUNGE  OF  ,I1,8H  EVECTOR , 2X , F7. 3 ) 

213 

66 

WRITE(6,215)J,TRE 

215 

215 

FORMAT ( 1 1 H  TREND  OF  ,I1,8H  EVECTOR , 2X * F8. 3 ) 

216 

WR I T£ ( 6, 2 10 ) V ( 1 , 1 ) 

217 

210 

FORMAT ( / 1 9H  LEAST  LATENT  ROOT , 2X , F 1 1 . 5 ) 

220 

T1=DSQRT (CHISQ/ (V(3,3)-V(l,l))) 

221 

T2=DSQRT (CHISQ/ (  V  (  2  »  2 )— V (1,1))  ) 

222 

T 1 T=T 1 /DSQRT ( l.D0-Tl**2) 

223 

THET1=DATAN( TIT )/CF 

224 

T2T  =  T2/DSQRT ( l.D0-T2**2) 

225 

226 


THET2=DATAN( T2T )/CF 
WRITE(6,212)THET1 


227 

212 

FORMAT { 22H  RADIUS 

OF 

CONFIDENCE 

230 

WRIT£(6,212)THET2 

231 

50 

CONTINUE 

232 

NMTRX=NMTRX+ 1 

233 

I F ( NMTRX.LT. MIX  )G0 

TO 

111 

236 

CALL  EXIT 

237 

END 

« 
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I  SN 

0 
1 
2 


34 

35 

36 

41 

42 

43 

44 
47 

50 

51 
52_ 

53 

54 

55 

56 

57 
60 


62 


63 


64 


65 

66 
67 

70 

71 
74 
77 

1102 
103 
li  0  4 


SOURCE  STATEMENT 
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3 

4 

5 

6 
7 

12 

13 

14 

15 

16 

23 

24 

25 

26 
27 

32 

33 


SIBFTC  JACOBI  DECK 

SUBROUTINE  JACOBI ( N , MATRI X , VECTOR , TOLERC , NORM ) 

_ DOUBLE,  PRECISION  SINE, COSINE, T S 1 , T S 2, LAMDA, MU, OMEGA , TOLERC , 


1 

2 


I  MATRIX (20, 20) , VECTOR ( 20 , 20 ) ,T$3 

N  =  3 

N0RM=2 

T0LERC=0. 000000000000000 
INTEGER  I, J,P,Q,N,NMl,PPlt DONE ,NORM 
_I£JL  T 0 L E R C . E 0,0. 00 00000000 00)  T0LERC=0.00 00000000001 
NM 1  =N— 1 
DO  1  1=1, N 
DO  1  J=  1 , N 

VECTOR ( I , J ) =0.000000000000000 

IF (  I  .EQ. J )  VECTOR ( I , I )  =  1. 000000000000000 

_D  0  N  E  =  0 
DO  5  P= 1 , NM 1 
PP1=P+1 
DO  5  Q=PP1,N 


IF ( DABS ( MATRI X ( P,Q)  )  .LE. TOLERC ) 
DQNE= 1 

LAMDA=-MATRIX(P,Q) 


GO  TO  5 


MU=.5OO0C0OO0OOOOOOO*{MATRIX(P,P)-MATRIX(Q,Q) ) 
OMEGA=LAMDA/DSQRT (LAMDA*LAMDA+MU*MU ) 

I F ( MU.LT. 0.0000 00000000000)  GMEGA=-OMEGA 
SIN E=OM EGA/ DSQRT(2. 00000000000000 0+2.  000000000000000* 
1  DSQRT (  1.000000000000000-OMEGA*OMEGA) ) 

__COS INE  =  DSORT( 1 . OOOOQQQQQOQOOOO-SI NE*SI NE  ) _ 

DO  4  1=1, N 

IF(  I.EQ.P.OR. I.EQ.Q)  GO  TO  3 
TS1=C0SINE*MATRIX(P, I ) -S I NE*MATR I X ( Q, I ) 
TS2=SINE*MATRIX(P,I ) +COS I NE*MATR IX ( Q, I ) 

MAT  R  I X ( P  » I )  =  TS1 

M ATR  I X ( Q , I )=TS2 _ 


MATRIX ( I,P)=TS1 
MATR  I  X ( I , Q ) =TS2 

TS 1  =  VECTOR ( I ,P) * COS  I NE- VECTOR (  I,Q)*SINE 
T  S2  =  VECT  OR ( I , P ) *S I NE  + VECTOR ( I ,Q)*COSIN£ 
VECTOR ( I ,P)=TS1 
VEC  TOR ( I , Q ) =TS2 


TS 1  =  MATR I X ( P , P )*COSINE» CO SINE- 2. 000000000000000* MATRIX(P,Q)* 
1  COS INE*SI NE  +  MATRIX(Q,Q)*SINE*SINE 

TS2  =  MATRIX( P,P) *SINE*SINE  +  2. 000000000000000 *MATR I  X ( P , Q ) * 

1  COS  I NE*S IN  E  +  MATRIX (Q,Q) *COS I NE*COS INE 

TS3=(MATRIX(P,P  i-MATRI  X  (  Q ,  Q )  )  *SINE*C0S1NE  MATRIX(P,Q)* 

X  (COSINE  *  C 0  SI NE-S INE*  S INE) _ 

MATRIX! P,P)=TS1 
MATRIX! Q,Q)=TS2 
MATRIX! P,Q)=TS3 
MATRIX! Q, P)=TS3 
CONTINUE 

IF! DONE. NE.O)  GO  TO  2 _ 

IF! NORM.NE. 1)  GO  TO  8 
DO  7  J= 1  *  N 
TS1=VECT0R( 1 , J) 

DO  6  1=2, N 


. 

. 


. 
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105 

6 

IF (DABS( TS1 ) .LT. DABS (VECTOR! I , J) ) ) T  SI  =  VECTOR(I,J) 

111 

DO  7  1  =  1,  N 

112 

7 

VECTOR ( I , J  )  = VE CTOR { I , J )/TSl 

115 

RETURN 

116 

8 

DO  11  J  =  1 ,  N 

117 

TS1  =  0.000000000000000 

120 

DO  9  1=1, N 

121 

9 

TS1  =T S  1  +  VECTORd  ,J)*VECTOR(  I, J) 

123 

TS1  =  l.COOOOOOOOOOOOOO/DSQRT (TSi) 

124 

DO  10  I  =  1 ,  N 

125 

10 

VECTOR ( I , J  )  =  VECTOR ( I , J ) *TS1 

127 

11 

CONTINUE 

131  RETURN 

132  END 


_ 74 


450008 

cCT  PROGRAM  IS  BEING  ENTERED 

IBLCR  —  JOB 

INTO  STORAGE. 

vc  tors 

IX  NO.  1 

[I CN  NUMBER  3  73 

OF  READINGS  IC 

V- 

i 

jrorn 

3 

— 

II  CN  NUMBER  3  26 
OF  READINGS  10 

r I  CN  NUMBER  32  7 

OF  READINGS  10 

- - - — — - - - - - - 

- - - 

OF  STATIONS  3 

• 

;  UNDERFLOW  AT 

*GE  OF  1  EVECTGR 

X  OF  1  EVECTGR 

020467 

12.650 

316.524 

— _ — - - . — - - 

,GE  OF  2  EVECTGR 
,0  OF  2  EVECTGR 

21.550 

231.609 

>T  LATENT  ROOT 

0. 37171 

US  OF  CGNFIDENCE 

1.286 

US  CP  CONFIDENCE  7.846 


75 


•  j.L 

negative  station  number  signifies  the  p+1  station  begins  the  calculation  of  the 

matrix. 

Use  of  the  program 

The  user  should  prepare  a  deck  with  the  following: 

(1)  The  standard  control  cards  for  the  IBSYS  monitor-- 'at  the  University  of  Alberta 
these  consist  of  a  j>  IBSYS  card,  a  JOB  card,  a  j^TIME  card  five  minutes  is  ample 
time  to  compute  ten  matrices,  the  number  of  lines  of  output  will  depend  on  the  number  of 
matrices  (2  minutes  is  sufficient  if  a  relocatable  binary  deck  is  used). 

(2)  A  source  deck  or  the  corresponding  relocatable  binary  deck. 

(3)  A  ^ ENTRY  card. 

(4)  The  input  data  as  written  up  under  restrictions. 

Check  Out 

The  program  has  been  used  to  divide  a  large  anticline  into  cylindrical  portions. 
About  20  such  portions  were  determined.  A  rough  check  (to  within  10°)  on  these 
results  was  made  graphically.  These  checks  tallied.  The  internal  consistency  of  the 
results  can  also  be  checked  roughly.  This  proved  satisfactory.  There  is  no  way  at  the 
moment  of  making  more  accurate  checks  on  these  results. 

The  sub  routine  (JACOBI)  to  compute  the  eigen  values  of  the  matrices  was  written 
by  G.  Jackson  of  the  Department  of  Computing  Science  in  Fortran  IV  (double  precision). 
(A  program  since  published  in  the  SHARE  library  SDA  3202  01,  Real  Symmetric  Matrix 
Eigen  System  Solver,  supersedes  Jackson's  program).  A  print  out  of  JACOBI  is  included 
with  the  program  to  solve  W. 

2)  To  calculate  the  axes  of  folds  assumed  to  be  cylindrical. 


Purpose 


This  program  is  a  particular  case  not  covered  by  the  earlier  programs — N.  is  1  in 


all  cases. 
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SOURCE  STATEMENT 


CRT RAN  SOURCE  LIST 


SN 

'  C 

1 


3 

t; 

7 


IGF 


T 


l  iv  A  I N  NOLIST  ,  NOD  ECK 
DOUBLE  PRECISION  V ( 20 , 20 ) , S ( 20 , 20 ) , D ( 3 ) , P I f RZ , AL , AM , AN , TR , CF , PL , 

1  A?RE  _»DLIMtSLIM,A,B*C*6»RfY,Zt TRE , APLU, PLUtTOLERC 

N  M  t  R  X  =  0  - - ~ 

Q  F  API  S  _  ')  1  \  i  i  T  V 


21 

1 1 1 


R  fc  A  C  (  5 ,  2  1  )  i  1 1  X 
FORMAT (IX, 13, IX) 
CG22J= 1 , 3 
C  C  2  2  I  =  1 ,  3 


110 

22 

V (  I , J ) -DC  0 

13 

lU 

13 

P  I  =3  . 141  5  9  2  6  535  8  979 
RZ= 1 .00/ 180.00 

C  F  =  P  I  *  R  Z 

lo  COMpl'fo.  & 
\e&olve  ffce. 

lo 

32 

i\  =  0 

17 

21 

4  1 

READ (5,41  ) NSTN 

FORMAT ( IX , 13 } 

• 

Matrix^  U 

22 

AL=CCC 

23 

A  Fi  =  C  C  0 

• 

24 

A  N  =  0  b  0 

23 

READ  (5  ,51  ) 

LLIM,NLIM 

30 

5  1 

FORMAT ( IX , 13, IX , 13) 

3 1 

'X  ^ 

i\  =  N  +  1 

32 

READ  (5  ,52  ) 

NSTR , ND I P , ND I R 

3  6 

52 

FORMAT ( IX, 13,  IX  ,  12,  IX 

,11) 

37 

IF(NSTR.EQ.999)G0  TO 

1 1 

42 

PL=90—  NDIP 

43 

I  F ( N  D I R— 1 ) 3 , 2 , 3 

• 

44 

2 

TR  =  NS  TR  +  9  0 

43 

GO  TO  5 

40 

3 

T  R  =  N'  S  T  R—  9  0 

47 

I F  (  T  R  )  4 , 5 , 5 

j 

50 

4 

TR=TR+36C .DO 

51 

5 

CLIM=LLIM 

52 

I  F { TR-DLIM )6, 10 , 10 

53 

6 

SL  IM=NLIN 

54 

IF(SLIM-CLIM)8,9,9 

v  D 

8 

IF ( TR-90.CG) 10, 10,9 

36 

9 

TR=TR+1RC.D0 

57 

P  L  =  P  L  *  C  F 

50 

TR=TR*CF 

51 

E=-CS I N ( TR ) 

j2 

GO  TO  101 

.3 

10 

PL=PL*CF 

54 

TR=TR*CF 

55 

E  =  DS  IN { TR ) 

:o 

101 

A  =  C C G S ( PL  ) 

J 

B  =  DS  IN ( PL ) 

7C 

C  =  DCGS (TR  ) 

n 

12 

V ( 1 , 1 )=V ( 1 , 1 ) +A*A*C*C 
V  (  1 , 2-)=V  (  1 , 2  )  +A#A*C*E 

—  — - —  -  — - -  — 

V ( 1 , 3 ) =V ( 1 t 3 )+A«C * B 
V  (  2 , 2  )-V  (  2T3T+X* E * 3 


■ 


. 


ISN 

10  i 
102 
iv 

12 

„  hOKTRAN  SOURCE  LIST  MAIN 

SuURCE  STATEMENT 

W  R  I  T  E  (  6  , 12) N  S TN 

FORMAT  ( //  15H  STATION  NUMB  ER  ,  9X  ,  I  3  ) 

'a  R  I  T  E  (  6 , 1  5  )  N 

K 

1  5 

EU  IT  a i  ( 1 8H  NO.  OF  RE AD  I NGS , 6X , I 3 ) 

105 

13 

V { 2, 1  )=V (  1 ,2  ) 

106 

V( 3, 1 )=V( 1 , 3) 

1C7 

V( 3t2)=V(2,3) 

110 

CALL  JACOBI (  N  ,  V  ,  S »  TOLERC , NORM ) 

ill 

L0310L=1,2 

1? 

C031CK-1 , 2  ' 

i  13 

I i  ( DABS ( V ( K  ,  K ) )  .  LE. OABS ( V { K+l ♦ K+l ) )  )GO  TO  310 

116 

Y  =  V  (  K  ,  K  ) 

117 

V(K,K )=V(K+1,K+1  ) 

20 

C  0  3  1 C  I  =  1 , 3 

1 2 1 

V(K+1,K+1 ) =  Y 

[22 

Z=S(K,  I  )  "  ■  . 

.23 

S ( K, I )  =  S (K+l,  I ) 

S { K+l ,  I  )  =  Z 

,23 

310 

CENT  I NUE 

.31 

.  C  G  3  1 1 M  =  1  ,  3 

32 

311 

C  (  M  )  =  S  {  1  t  M  ) 

34 

IF(C(1))36,40,36 

35 

36 

ATR  E=D ( 2)/D( 1) 

36 

TRE  =  CAT  AN ( ATRE) /CF 

37 

GO  TO  37 

40 

40 

IF(C(2))43,42,42 

41 

43 

TRE=9C . DC 

42 

GO  TU  42C 

43 

42 

TRE=27O.C0 

44 

420 

A  P  L  U  =  0 ( 3 ) /DSORT ( 1 .DO-D ( 3 ) **2 ) 

45 

P  L  U  =  C  A  T AN ( APLU ) 

46 

P L U  =  D A B S ( PLU/CF  ) 

47 

GO  TO  65 

50 

37 

APLU  =  C( 3)  /DSQRT (l.DG-C(3)**2) 

31 

PLU  =  CAT  AN ( APLU) 

52 

PLU=DABS ( PLU/CF ) 

53 

IF ( C ( 3  )  )  6  i  ,60,60 

54 

60 

C ( 1  )  =-D ( 1 ) 

53 

C (2  )=-D( 2) 

36 

61 

I  F ( D ( 2 ) ) 63  f  62  t  62 

57 

62 

IF ( D ( 1 ) ) 64 , 65,65 

60 

6  3 

TRE= 180 . CO+TRE 

61 

GO  TO  65 

6  2 

64 

TRE=360. CO+TRE 

63 

63 

WRITE(6,214)PLU 

34 

214 

FORMAT ( 2 1 H  PLUNGE  OF  FOLD  AX  I S , 2X , F7 . 3 ) 

k  r: 

aR ITc ( 6, 215 ) TRE 

6 

215 

FORMAT  (  20H-  TREND  OF  FOLD  AX  I  S  , 2X  ,  F 8 . 3  ) 

'  7 

aR ITE  ( 6, 2  10  )  V ( 1 ,  1  ) 

/o 

210 

FORMAT (  1 9H  LEAST  LATENT  ROO T , 2X , F 9 . 3 ) 

71 

i\MTRX  =  NMTRX  +  1 

72 

IF ( NMTRX .LT.MIX )GO  TO  111 

75 

50 

CONTINUE 

7 

CALL  EXIT 

/ , 

END 

» 


' 
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Method 


i  i  i 


Watson  (I960)  leads  to  the  conclusion  that  the  fold  axis  of  the  readings  (I.,  m.,  n.) 
is  the  eigen  vector  associated  with  the  least  eigen  value  of  the  matrix  U, 

JTl,  ^El{n, 

\;M-  2T  Yv\*~  *2^  fni  h, 

.>  1,-n,  KTh,  2  h'A 


u  = 


Watson  (1965)  suggested  that  the  least  eigen  value  might  be  the  basis  of  a  test  of 
coplanarity, 

(least  eigen  value  of  U) 

when  N  is  the  total  number  of  readings. 

Restrictions 

The  input  data  must  be  in  the  form  suggested  in  Program  1  with  the  proviso 
that  "station"  now  be  interpreted  as  the  population  of  N  readings  whose  axis  is  to  be 
determined. 

Negative  station  numbers  need  not  be  used. 

Use  of  the  program 

As  in  program  1 . 

3)  Calculation  of  the  Coplanarity  Statistic 

Purpose 

The  purpose  of  this  program  is  to  calculate  a  coplanarity  statistic  to  test  the  nul 


hypothesis  that  a  number  of  samples  of  planar  structures  or  lineations  come  from  a  single 
population. 
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SOURCE  STATEMENT 


bo _ : _ _ _ _ 

FORTRAN  SOURCE  LIST 


|CSN 


0  $  I BFTC  MFPL  NOLIST, DECK 

1  PI  =  3. 1415926 

2  R  180  =  1 ./ 180. 


3 

4 

5 

7 

10 

11 

53 

13 

CF=PI*R180 

NRT  =  0 

READ ( 5 , 53 ) NR 

FORMAT ( IX, 13, IX ) 

NRT =NRT+ 1 

T  STN-0 

12 

13 

14 
13 
16 

17 

68 

T  K  V  =  0  — j—  |  |  L 

nt=o  lo  calculable 

* SG=0  tu  Copla  v\c\  ki  fv 

ts  a=o  Si  aTisr \c ,  C.  •  S  . 

READ!5,53)NSTN 

21 

30 

FORMAT! IX, 13) 

22 

IFINSTN.LT. 0)  GO  TO  69 

25 

WRI TE! 6,60  )NSTN 

26 

60 

FORMAT ( /8H  STATI0NI4) 

27 

TSTN-T  STN+1* 

30 

READ! 5,51)LLIM,NLIM 

33 

51 

FORMAT! IX, 13, IX, 13) 

34 

SA=0.0 

35 

S8=0 . 0 

36 

SG  =  0 . 0 

37 

NCDS=0 

40 

1 

READ(5,52)NSTRf NDIP,NDIR 

44 

52 

FORMAT!  IX, 13,  IX, 12, IX, ID 

45 

IF(NSTR.LT.O)  GO  TO  14 

50 

NCDS=NCDS+1 

51 

PL=90-NDIP 

52 

IF(NDIR-1)3,2,3 

53 

2 

TR=NSTR+90 

54 

GO  TO  5 

, 

55 

3 

TR=NSTR-90 

56 

IF! TR)4,5,5 

57 

4 

TR=TR+360. 

50 

5 

DL I M=LL  IM 

51 

IF(TR-DLIM)6,10»10 

52 

6 

SL I M  =  NL  I M 

, 

53 

IF! SLIM—  DLIM)8,9,9 

54 

8 

IF(TR— 90. )10,10,9 

55 

9 

TR=TR* 180. 

56 

PL=PL*CF 

* 

57 

TR=TR*CF 

JO 

A=COS!TR)*COS!PL) 

n 

B=SIN(TR)*COS(PL) 

'2 

G=-SIN! PL) 

r3 

GO  TO  11 

r4 

10 

PL=  PL*CF 

r5 

TR=TR*CF 

76 

A=C0S ! TR ) *C0S ! P  L ) 

n 

B  =  S  IN ( TR)*COS(PL) 

G=SIN! PL) 

• . 

11 

SA=SA+A 

.  . . 


'  ■ 


_  81 

FORTRAN  SOURCE  LIST  MFPL 

SOURCE  STATEMENT 

SB=SB+B 
12  SG=SG+G 

_ GO  TO  1  _ 

14  R=SQRT { ( SA*S A ) +  (SB* SB )  +  ( SG* S  GT ) 

WR I TE ( 6  »  65 ) NCOS 

65  FORMAT ( / 13H  MEASUREMENTSI4) 

WRITE(6,59)R 

59  FORMAT ( 14H  VECTOR  LENGTH, 3X,F9.4) 

T  S  A  =  TS  A-a-S  A 

TSB=T  SB+SB  — 

TSG  =  T  SG+SG 
NT  =NT+NCDS 
TRV=TRV+R 
GO  TO  68 

69  RV  =  SQRT ( { TS A*TS A )  +  ( TS8*TSB )  ♦  (  T SG*TSG ) ) 

WR  I  T  E  (  6  9  6  1 )  R  V 

61  FORMAT ( / / 17H  RESULTANT  VECTOR , F9.4 ) 

FNT=NT 

F=(FNT-TSTN)MTRV-RV)/( ( TSTN-1 . ) * ( FNT-TRV ) ) 

WRIT6C6,62)F 

62  FORMAK/lOH  STATIST  IC,7X,F9. 4) 

IF  <  N  R T . L  T . NR) GO  TO  1 3 

CALL  EXIT 
END 
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4  6  C  0  0  5 


IblCR  —  JOB  MFPL 


JECT  PROGRAM  IS  BEING  ENTERED  INTO  STORAGE. 
I G f \  373 


SI  R  E  N  ENTS  10 
CR  LENGTH  9.9828 

TICK  32  6 


ilREMENTS  10 

tcr  length;  9.9728 


I C  N  32  7 


^  fa  '  1 

F/J  b  - 


row\ 


REF,  ENTS  10 
I R  LENGTH 


9.9738 


LTANT  VECTOR  29.4640 


{STIC  _ 

re  i-r 


91 . 8  76  6 

'ca  v\  l~ 


✓  ence-  ^cov*\  Co  ✓  i  lCj 


83 


Method 

The  lineations  or  normals  to  the  planar  structures  are  assumed  to  be  distributed 
about  the  estimated  means  of  their  sample  according  to  the  probability  density  function 

P(0), 

P  (9)  =  B  e  COS^  (Fisher,  1953) 

Watson  (1956)  has  shown  that  the  statistic 

2>;-F  .  ZR,  -d 

(P-V  Z*;  •'  ' 

referred  to  F-ratio  tables  with  2(p~l)  and  2(  N.“p)  degrees  of  freedom  is  a  test  of  the 

t*h 

above  hypothesis.  The  i  sample  contains  N.  readings  whose  vector  sum  has  a  length 
R..  R  is  the  length  of  the  vector  sum  of  the  resultants  of  the  readings  and  p  is  the 
number  of  samples. 

The  program  outputs  values  of  R.,  N  and  the  coplanarity  statistic. 

Restrictions 

The  input  data  for  this  program  must  be  in  the  following  form: 

Card  1 

The  number  of  separate  tests  of  coplanarity  to  be  made  in  this  run  is 
punched  in  2,  3,  4. 

Card  2 

A  3  digit  number  identifying  the  sample  to  the  programmer  is  punched  in 

2,  3,  4. 

Card  3  and  the  data  cards  are  as  described  in  the  programs  for  calculating  the 
generatrix  of  folds. 

The  last  reading  within  a  sample  is  followed  by  a  data  card  with  a  negative 
trend  punched  in  it  (in  the  format  of  Card  3). 
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The  last  reading  within  the  group  of  samples  whose  coplanarity  is  being  tested 
is  followed  by  a  negative  trend  and  a  negative  identifying  code  (in  the  format  of  Card  2). 
Use  of  the  program 


The  user  of  the  program  should  prepare  a  deck  with  the  following: 

(l)  The  standard  control  cards  for  the^lBSYS  monitor-—at  the  University  of 


Alberta  these  consist  of  an  IBSYS  card,  a  ^jLlOB  card  and  a  ^TIME  card — -thirty  seconds 


is  sufficient  for  the  usual  size  test. 

(2)  A  source  deck  of  the  corresponding  relocatable  binary  deck. 

(3)  A  jjENTRY  card. 

(4)  The  input  data  as  written  up  under  restrictions. 


Check  out 

The  values  of  R  and  R.  were  checked  by  Muecke's  mean  S  plane  program 
(Muecke,  1965).  A  calculation  of  the  statistic  from  the  R  values  by  hand  corresponded 
with  the  figures  produced  by  the  program. 

4a)  Preparation  of  Pole  Density  Diagrams 

The  program  given  below  is  a  translation  into  Fortran  IV  and  slight  modification 
(by  V.  Yanda  and  D.  Cruden)  of  the  program  written  in  Fortran  II  by  V.  Yanda  and 
G.  Muecke  (Muecke,  1965,  pp.  18-30).  The  program  is  discussed  exhaustively  there 
and  only  modifications  will  be  discussed  here. 

Parameter  Cards 
Card  1 

A  1  digit  number  is  punched  in  column  3.  This  determines  the  size  of  the  cross 
section  of  the  counting  cone  (1%  -  9%  of  the  area  of  the  hemisphere)  which  surrounds 
the  grid  points. 
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Card  2 

A  3  digit  number  which  identifies  the  station  to  the  programmer  is  punched  in 


If  the  data  is  in  the  form  of  the  plunges  and  trends  of  lineations,  1,  is  punched  in 
column  6. 


A  negative  station  number  indicates  the  end  of  a  run  of  several  stations. 

Card  3  -  data  cards 

Data  cards  are  in  Muecke's  original  format  except  that  a  negative  strike  now 
indicates  the  end  of  a  set  of  data  from  a  particular  station.  It  should  be  followed  by 
a  card  of  type  2. 

Use  of  the  program 

The  user  should  prepare  a  deck  with  the  following: 

(1)  The  standard  controle  cards  for  thejjBSYS  monitor— at  the  University  of  Alberta 
these  consist  of  an 

ient  to  process  300  readings  in  a  typical  station  if  a  relocatable  binary  deck  is  used, 
each  station  requires  120  lines  of  output. 

(2)  A  source  deck  or  relocatable  binar  deck. 

(3)  A  ^EN TRY  card. 

(4)  The  point  density  counter  given  by  Muecke  (1965,  Appendix  1). 

(5)  The  input  data  with  the  appropriate  parameter  cards. 
b)  This  program  for  the  preparation  of  pole  density  diagrams  has  been  modified  in  a  way 

described  in  Appendix  5.. 

Parameter  cards 
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FORTRAN  SOURCE  LIST 


1SN  SOURCE  STATEMENT 


0 

1 

o 

L. 

3 

4 
6 

11 

13 

14 

15 

16 
17 
21 
22 

23 

24 
26 
26 

30 

31 

32 

35 

40 

41 

42 
46 

51 

52 

53 

56 

57 
60 
61 
62 
65 

70 

71 

72 

73 

74 

75 

76 

77 
00 
01 
02 
03 
04 
03 
06 
07 

12 
13 

i  r 


i  I  8FTC  CCOUNT  NGOECK , NOL I  ST 

Cl  MENS  I  ON  COSA(  33  3)  , C0S6(  3  33  ) , COSG(  33  3) , NO (333) 
DIMENSION  A  (  2000  )  ,B  (  2000  )  ,G(  20C0  ) 

CO  2  1=1,333 

RE AC ( 5, 51  )COSA(  I  ) ,COSB ( I ) ,COSG (I) , J 
I F (  I  .NE. J  )  GO  TO  23 
2  CONTINUE 


3 

4 


301 

300 


P 1=3. 1415926 
R180=l./180. 
CF=PI*R180 


CO  4  1=1,333 
ND(  I  )  =  0 
NCDS=0 
NZ  =  0 
TNO~0 
I  X  =  1 


R£AD(5,3C0)QTA,QTB,QTC,0TD 
READ ( 5 , 3C1 ) THA, THB, THC ,THD,K,DCF 
FORMAT ( 1X,4F5.0 ,2X, 14, F5.2) 
FORMAT ( IX, 4F  11.8} 

REAC ( 5, 5C )NSTN, NZ 
IFtNSTN.LT.O  )CALL  EXIT 


NZ=NZ+ 1 


GO  TO  (5,10)  , NZ 

5  READ  (5  ,53  )  NSTR , NOI P , NO  I R 

IF(NSTR.LT.O)GO  TO  12 

_ NC_DS_=NCDS+1  _ _ _ 

PL  =  90—ND  IP 


IF (NDIR.NE.  1)  GO  TO  8 
TR=NSTR+90 
C-C  TO  11 
8  TR=NSTR-9C 

_ GO  TO  11  _ 

10  READ  (5,535)  NTR ,NPL 
I F ( NTR.LT .0 ) GO  TO  12 
NCDS=NCDS+ 1 
TR=NTR 
PL=NPL 


11  PL=PL*CF 
TR=TR*CF 


I  =  N  C  C  S 

A (  I  )=COS  ( TR ) *COS  (PL) 

B (  I  )  =  S  I N  ( TR ) *COS  (PL) 

G (  I )  =  S I N  (PL) 

GO  TO  (  5  ,  1 0  )  ,  N Z _ _ _ 

12  CONTINUE 

FNCDS=NCCS 
RNCDS=100./FNCDS 
p  M  A  X  =  0  . 

CO  124  1=1, NCOS 

IF ( PMAX . GE . G ( I ) )  GO  TO  124 

P  M  A  X  =  G  (  I  ) 

124  CONTINUE 

PMAX=PMAX+QTA 
DO  15  1=1,333 
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17 

CSG=COSG ( I ) 

20 

C  E  N  S  =  0  . 

n 

I  F ( PMAX.LE .CSG) 

GO  TO  15 

..4 

CSA=COSA ( I ) 

2b 

CS8=C0SB ( I ) 

26 

CO  14  J=  1 , NCDS 

27 

GEE=G ( J ) 

30 

I  F { ABS ( GEE-CSG) • 

GT. QTA ) GO  TO  14 

33 

R  E  X  =  A  8  S  (  C  S  A  *  A  (  J )+CSB*B(J)+CSG*GEE) 

34 

IF ( QTA.GT.REX )GQ 

TO  14 

37 

CENS  =  CEN  S  +  TFIA 

40 

IF(QTB.GT.REX)GO 

TO  14 

43 

DENS  =  DtNS+  TUB 

44 

I F ( QTC .GT.REX )GO 

TO  14 

47 

DENS=CENS+THC 

50 

IF ( QTD. GT.REX ) GO 

TO  14 

53 

CENS=DENS+THD 

54 

14 

CCNTINUE 

56 

13 

NO (  I  )  =  0 E N S 

60 

KL  I  K= 1 

61 

no  25  1=1,333 

62 

V  N  C  =  N  D  (  I  ) 

63 

25 

TNC=TND+VND 

65 

IF ( TND.LT.20. 

«F NCOS ) GO  TO  302 

70 

TND=FNCDS*20. 

71 

302 

CONTINUE 

72 

20 

CO  21  1=1,333 

73 

X  N  0  =  N  D  (  I  ) 

74 

21 

ND (  I  )  =XND*  100 

./TND+.5 

76 

16 

WRITE  (6 

♦  5  8 

)  N  D  (  1 ) 

77 

WRITE  (6 

,66 

)  ND( 318) ,ND(319) ,ND(320) ,ND(321) 

!0 

WRITE  (6 

,68 

)  NSTN 

)1 

WRITE  (6 

,59 

)  ( ND ( I ) , 1=2 , 10 ) 

16 

WRITE  (6 

,57 

) 

)7 

WRITE  (6 

,57 

) 

0 

WRITE  (6 

,60 

)  ( ND ( I ) ,1=11,23) 

5 

WRITE  (6 

,57 

) 

6 

WRITE  (6 

,57 

) 

7 

WRITE  (6 

,61 

)  (ND( I ), 1=24,38) 

4 

WRITE  (6 

,57 

) 

5 

WRITE  (6 

,57 

) 

6 

WRITE  (6 

,62 

)  (ND (I), 1=39 ,55) 

3 

WRITE  (6 

,57 

) 

4 

WRITE  (6 

,57 

) 

5 

WRITE  (6 

,62 

)  (ND( I ) ,1=56,72) 

2 

WRITE  (6 

,57 

) 

3 

WRITE  (6 

,57 

) 

4 

WRITE  (6 

,  6  3 

)  (ND (I), I =73*91) 

1 

WRITE  (6 

,57 

) 

2 

WRITE  (6 

,57 

) 

3 

WRITE  (6 

,63 

)  (ND( I ) , 1=92, 110) 

0 

WRITE  (6 

,57 

) 

1 

WRITE  (6 

,57 

) 

~> 

WRITE  (6 

,63 

)  (ND( I ) , 1=111, 129) 

. 

WRITE  (6 

,67 

)  ND ( 322 ) , ND ( 32  3 ) 

. 


— ■  % 


- 


• 

. 

* 

94 _ 

FORTRAN  SOURCE  LIST  DCOUNT 


SN 

SOURCE 

STATEMENT 

70 

WRITE 

(  6 

,57 

) 

71 

WRITE 

(  6 

,63 

)  (ND( I ) , 1=130, 148) 

-'b 

WR  I  TE 

(6 

,67 

)  ND( 324) ,ND(325) 

.7 

W  RITE 

(  6 

,57 

) 

00 

WRITE 

(  6 

,64 

)  (ND ( I ) , I =149, 169) 

05 

WRITE 

(  6 

,57 

) 

06 

WRITE 

(  6 

,67 

)  ND( 326) ,ND( 327) 

07 

WR  I  TE 

(6 

,63 

)  (ND( I ) , 1=170, 188) 

14 

WRITE 

(6 

,57 

) 

15 

WRITE 

(  6 

,67 

)  ND( 328) , NO (329) 

16 

W  R  I  T  E 

(6 

,63 

)  (ND( I ), 1=189,207) 

23 

WRITE 

(  6 

,57 

) 

24 

WRITE 

(  6 

,57 

) 

25 

WRITE 

(6 

,63 

)  ( ND ( I ) ,1=208,226) 

32 

W  RITE 

(  6 

,57 

) 

33 

WRITE 

(6 

,57 

) 

i 

34 

WRITE 

(  6 

,63 

)  ( ND ( I ) ,1=227,245) 

41 

WRITE 

(  6 

,57 

) 

42 

WRITE 

(  6 

,57 

) 

43 

WRITE 

(  6 

,62 

)  ( ND ( I ) , 1=246,262) 

30 

WRITE 

(6 

,57 

) 

51 

W  R  I  T  E 

(  6 

,57 

) 

WRITE 

(  6 

,62 

)  ( ND ( I ) ,1=263,279) 

>7 

WRITE 

(  6 

,57 

) 

50 

WRITE 

(  6 

,57 

) 

>1 

WRITE 

(  6 

,61 

)  ( ND ( I ) ,1=280,294) 

)6 

WRITE 

(6 

,57 

) 

j1 

WRITE 

(  6 

,57 

) 

JO 

WRITE 

(  6 

,60 

)  ( ND ( I ) ,1=295,307) 

f5 

WRITE (6, 54) NCOS 

7b 

WRITE(6,69)K 

7 

WR  I  TE 

(  6 

,59 

)  ( ND { I ) ,1=308,316) 

!4 

WRITE 

(6 

,57 

) 

15 

WRITE 

(6 

,66 

)  ND ( 330 ) , ND ( 33 1 ) ,ND(332) ,ND(333) 

16 

WRITE 

(  6 

,65 

)  ND ( 3 17 ) 

7 

GO  TO 

( 3C3  »  22 ) , KL IK 

0 

303 

C03G4  I 

=1,233 

1 

XND  =  ND (  I  ) 

2 

304 

R  D  (  I  )  = 

XND*FNCDS*DCF/200.+.5 

4  K  L I K  =  2 

5  GO  TO  16 

6  22  GO  TU  3 

7  23  WRITE (6, 71) 

0  CALL  EXIT 

1  24  WRITE (6, 72) _ 

2  CALL  EXIT 

3  5C  FORMAT ( 1X» I  3, IX  > ID 

4  31  F0RMAT(.7Xf  3F11. 8,  12X,  13) 

5  32  FORMAT (13) 

53  •  FORMAT  C*1X  ,  I  3,  IX  ,  12,  IX,  1 1 ,  IX,  15  ) 

1  335  FORMAT  f  IX, I  3,  IX  ,  12, 3X,  15) _ ; _ 

0  57  FORMAT ( IX  ) 

1  53  FORMAT ( 21H1DENSITY  D  I  STR I  BUT  I  ON , 29X , I  5 ) 

0  59  FORMAT ( 3CX , 9 1  5 ) 

->  60  FORMAT  (  2CX,  1315  ) 


' 
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. 
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61  FORMAT (  1  5X  T  1516  ) 

62  FORMAT (  1CX  ,  1715  ) 

63  FORMAT ( 5Xf 1915) 

64  FORMAT  (21 15) 

69  FORMAT ( 4H  K=,I4) 

5A  FORMAT ( 2X, 14, 9H  READINGS) 

65  FORMAT ( 5CX , 15) 

c 6  FORMAT { 39X , I  5, 1  6, 1 10 , 16) 

67  FORMAT (  I6,93X,I5) 

6o  FORMAT ( 1 2 h  STATION  NO,, 14) 

70  FORMAT ( 4 H  PER,I3,13H  PERCENT  AREA) 

71  FORMAT ( 33H  ERROR=  CONSTANTS  OUT  OF  SEQUENCE) 

72  FORMAT (52H  ERROR=  PERCENTAGE  OF  AREA  CONSTANT  MISSING  OR  ZERO.) 

END 
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a  reading  af  that  grid  point  changes  value  are  punched  on  this  card.  They  start  at  the 
largest  angle  and  specify  its  direction  cosines  to  eight  figures.  Values  are  taken  from 


Table  5:1 . 


Card  2 


The  increments  of  each  weighting  are  punched  in  5,  10,  15  and  20.  The  value 
of  the  precision  parameter  assumed  in  these  calculations  is  punched  in  23  -  26  (no 


l 


ched  in  27  30 


decimal  figures).  -  w  ,  for  this  value  of  the  precision  parameter  is  pun 

/  . 

(to  2  decimal  figures). 

Cards  3  and  4  are  as  cards  2  and  3  of  the  unmodified  program. 

The  use  of  the  program  is  the  same  as  the  unmodified  program. 

Checkout 

The  program  was  checked  out  by  comparison  with  the  unmodified  program. 

5)  Determination  of  Mean  S  planes 

This  program  was  originally  written  by  Muecke  in  Fortran  II  and  a  full  acount  of 
it  is  contained  in  Muecke  (1965,  pp.  31-40).  The  program  was  translated  directly 
into  Fortran  IV  by  V.  Yanda  and  D.  Cruden.  As  it  stands,  the  program  is  rather 
inelegant  though  quite  functional. 

Method 


Fisher  (1953,  p.  296)  showed  that  the  best  estimate  (/A  ,JA  ,  "yO  of  the  true 
mean  direction  is  the  vector  sum  of  the  individual  directions  (  ^  .,  JA  'y'.).  That 


where 


>  - 


R 


-  (Z\,f  i  r(iRr 


The  aximuth  and  the  inclination  of  the  resultant  vector  are  determined  using  the  relation¬ 


ships. 
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where  D  is  the  declination  and  I  the  inclination  of  the  resultant  vector. 

The  accuracy  of  the  calculated  mean  direction  (\,M,  V)  can  be  estimated  by 


h  T)  -  ^  2^1  '  s  I  h  _l  =  I 

.  ^  X;  -IT 


the  angle  between  this  mean  and  the  true  mean.  Fisher  (1953,  p.  303)  showed  that, 
at  a  probability  level  of  (l-R),  the  true  mean  direction  of  the  population  lies  within  a 


circular  cone  about  the  resultant  vector  of  semi-vertical  angle 


where  N  is  the  total  number  of  vectors  included  in  the  analysis  and  the  precision  para¬ 


meter,  K,  is  greater  than  3.  If  K  is  greater  than  3, 


~n^“R 


Use  of  the  program 


The  user  should  prepare  a  deck  with  the  following: 


(1)  The  standard  control  cards  of  the  IBSYS  monitor— at  the  University  of  Alberta 


these  consist  of  a  IBSYS  card,  a  JOB  card,  a  TIME  card — two  minutes  is  ample  time 
for  ten  samples  of  ten  readings  on  the  IBM7040  (if  a  relocatable  binary  deck  is  used). 


lines  of  output  are  shown  in  a  sample  of  the  output. 

(2)  A  source  deck  or  corresponding  relocatable  binary  deck. 

(3)  A  ^ENTRY  card. 

(4)  Muecke's  parameter  cards  1  and  2  before  each  5. 

(5)  A  data  set. 

Check  out 

The  program  rewritten  in  Fortran  IV  was  checked  using  Muecke's  test  data.  The 
results  were  identical. 
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Jositnabi  ©  aw  et 


450''0?  ,V1  F  A  N  fracturp  PLANF 
5  ,  2  0  o 


DR  GO  ,  NOOECK 

|r  m^pl 

—  n  I  '-‘FNS  T  ON— A  (  500  )  ,  °.-(  500)  >  G  {  5-00  ) 
PT  =**.  ■»  AT  5  0  26 


(  MOD  IF  I  Er>) 


R) po  =  i ./IPO. 

^Prp  T*Pl  °  0 

RcrAn  (  *  *  R  0  )  N5TN,,Nr>FN,Nrnc  ,np 

TF (N5TN.LT. 0)  GO  TO  6  0 

p_  p  =  N  p _ _ _ 

Rc  ah  (  k  ,51  )  L  L  I  M  ,  N  L  I  M 

5  A  =  0 . 0 


o  r  =  n .  o 
c  G  =  0  .  C 

no  1 2  !  =  1 > NODS 

I-R-F-AA — (  - ,-5  ? - ;-N^TR-^Nn  I  P  ,N0  I  R 

Pi_  =  o^-.\jni  p 
! F ( NOT  P-i  ) 2,2  ,R 
TP  =  N  c  TP-fO  0 


GG  TO  c; 


TR  =m 5TR-Q  0 

— T-c--{— PR- )  A-,-5-,  5 - 

TP=TR+2A0. 

I  r>LIM  =  LL.IM 
I F ( TR-HLIM ) 6,10,10 
5  L  T  M  =  N  L  T  M 
I F { 5LIM-DL IM ) 9  ,  Q  ,  Q 
r“i~F-(-T-R — 0-0 0  ,-]  — 

T  R  =  T  P-r  1  «C  . 


PL  = 
7R  = 

A  (  T 
Q  (  t 

£4-4- 
GO 
PL  = 
TR  = 
A  (  I 
°  (  I 
G4-I- 


P  l  *  C  F 
7P*r  - 

)=r 05  (TR)*r05 
)=5TN  (TR)*rOS 

.).a^SLT44 — (-P-L-) - 

1  1 


TG 

I 


-  r 


( TR ) *C0S 

( TR ) *ros 


(PL  ) 
(PL  ) 


(PL  ) 
(  PL  ) 


1)e.tev  Minatioi'l 
-of-—  M&4  h — S— 

Plants 


ci=fA  +  4  (  I  ) 
co  =  CP  +  P  (  T  ) 

c  G  =  GG+O ( I  ) 

R= C0DT  ( ( S A*5A ) + ( 5P*5« ) + ( SG*5G) ) 
c  t  n  1 1  nj  _  c  r  /  R 

-T-A  fyC-— C  D  /-  A-A - — - — - 

R  F  T  N = A  7  A N  (GININ/SGRT  ( 1 . - 5 1 N I N*S I N I N )  ) 
RE  IN=PFIN/CF 
R  F  T  N  =  A  P  c  (  R  F  I  N  ) 

P  F  n  F  =  A  T  A  N  (TAOF) 

REHE=PEr>F/CF 

-IF  (  GG)  i-2  1  ,  122 ,1-2  2 - 

ca=- ( G  A ) 

C  p  __  (  C,D  ) 

I  F  (  c  A  )  1  c  ,  1  2  ,  ]  ? 

!  F  (  ^  )  i  5 , 1  7, 1  ~f 
o r n  ~ = i oo.+Rfnr 
GO.  _T0-  1  7 
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P^nF=?AO.+RFDP' 

rncrvrnc 

X  <  =  {  ros-i .  )  /  (  GDS-R  ) 

;TRJTf  (6  ,60  )  NSTN  »  NDFN 

WRITP  (A  ,6S  )  NCOS 

j—  ’  t  I  :  r-  (  6 —  — 6  “3- - ) — REDE - 

U/R IT-  (a  ,6?  )  RFI N 

’*■'  ^  T  T  F  (A  ,61  )  XK 

pr  o r= ( i ,-P  j 

P_  P^  _  A  "I 

0=1  ./  (  ^c-1  .  ) 

_od  -  ( r pc_p ) /p - 

PP= ( ( i . /P ) **Q ) -1  . 

9  R  =  RR*PP 

rOSAL=l .-PR 

F U i\  =  A R S  (  (  1  .  -fOSAL*rOSAL  )  /  (  fOSAL  )  ) 
FUN=SORT  (FUN) 

—A  L-=  A  Tan  -(-FUN- ) - - - — - 

A  L  =  A  L  /  r  F 

ncLTA=FQRT  ( ?.*( rnF-R ) /fDS ) 
rT A  =  O  F  L  T A/FQPT  (OOF) 

AFLTA  =  nci_T  A/fr 
r T  A  =  rT  A  f  C F 


-WO-J-T-P- 

-(A 

,-6-4 

)  P  E-R-r  ,  AL 

WR  I  Tc 

(  6 

,66 

)  DFL  T  A 

up  j  TF 

(  6 

,67 

)  FT  A 

GO  TO 

68 

0  A  L  L  F  X  I  T 

FORMAT ( IX, IF, IX, I?, IX, 13, IV, II  ) 

— F~0  R  M  A-T-pi- X-,- 1 -F-  ,-1-X-,-I  F-) - - : - 

FORMAT (  IX,  I  F  ,  1  x.  I?,l  X,  T1  ) 

FORMAT  (  //oh  fja.t  i  ONI  4 , 1  4H  M^AN  F  R  A  f  T  u  R  F  I  F  ) 

FORMAT  (?0H  PRFr'jc,jON  PAPAMCTFP  ,  1  3X  ,  FI  o.  F  ) 

F  0  P  M  AT  (  p  7  H  PFSI'LTANT  j  NO  L  T  N  A  T  I  ON  ,  1  OX  ,  F«  •  F  ) 

FORMAT (??H  Rf CULT  A NT  HFTL INAT TON , lox * FP . F ) 
if-0 j9-.vva-T-H-X-,-F-4  . 0-T4-X-,T?-F-H-t^-R-40&VF-FQAj FI  BENCE  R  A D I U 5-  F p  •  F -) 
FORMAT ( // I FH  MEASUREMENTS T?F ) 

FORMAT (14H  RMS  DEVIATION* 1 8X *FS. F ) 

FORMAT (1  PH  STANDARD  DEVI  AT  ION , 13X *F8 . 3 ) 

^ A  LL  FXIT 
FND 

l - M  F-P-L- 


• 

• 
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FI  /. 


t  J 


I  BLDR 


JOB  MFPL 


JcCT  PROGRAM  IS  BEING  ENTERED  INTO  STORAGE 


ION  37  3  N  E AN  FRACTURE  2 


oi'\C;"Li\  i  b 

LTANT  DECLINATION 
L  T  A \ T  INCLINATION 
ISICN  PARAMETER" 


10 

207.932 

_ 66.947 

321. 763 

PER  CENT  CONFIDENCE  RADIUS  2.116 
DEVIATION  3.365 

CARD  DEVIATION  1.064 


ION  326  MEAN  FRACTURt 


c  I 


U REN ENTS 

LTANT  DECLINATION 


10 

218.248 

_ 50.365  __ 

ISICN  PARA M E T E R  330 .  720 

PER  CENT  CONFIDENCE  RADIUS  2.659 


TANT  INCLINATION 


I Results 
jo*  stations 
Fi^,  3. 


DEVIATION 
CARD  DEVIATION 


4.227 

1.337 


I  ON  327  MEAN  FRACTURE  2 


UR  EM ENTS 

LTANT  DECLINATION 
.TANT  INCLINATION 


ISICN  PARAMETER  371.885 

PER  CENT  CONFIDENCE  RADIUS  2.507 
DEVIATION  3.986 

I.ARD  DEVIATION  1.261 

I F-D AT  A  ENCOUNTERED  ON  SYSTEM  INPUT  FILE. 


10 

201.941 

73.379 


I 


' 


i 
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APPENDIX  5 

THE  CONTOURING  OF  PI  DIAGRAMS 
The  relative  concentrations  of  readings  in  different  parts  of  a  pi  diagram  are 

often  best  displayed  by  contours.  These  concentrations  have  been  used,  Friedmann 

(1964,  pp.  468-469),  to  make  statistical  inferences  about  the  population  of  which 

the  readings  are  a  sample.  This  use  makes  it  clear  that  contouring  attempts  to 

estimate  the  probabilities  of  finding  readings  within  different  small  areas  on  the 

sphere.  It  is  implied  that  one  is  five  times  more  likely  to  find  a  point  within  the 

5%  contour  than  within  the  same  area  within  the  1%  contour. 

The  distribution  of  probability  is  customarily  investigated  by  Monte  Carlo 
methods  -  throwing  grids  on  circular  counters  across  the  equal  area  projection  on 
which  the  readings  have  been  plotted  and  counting  the  number  of  points  within 
each  cell.  The  number  of  readings  within  the  cell  is  then  an  estimate  of  the  prob¬ 
ability  of  finding  a  reading  within  a  small  area  at  the  centre  of  the  cell.  These 
probability  values  can  be  contoured  (usually  as  percentages  of  the  total  numbers  of 
points  plotted)  and,  it  is  argued,  the  distribution  of  probability  displayed. 

One  source  of  inaccuracy  in  this  method  has  been  pointed  out  by  Strand  (1944) 
Although  the  points  are  plotted  on  an  equal-area  net,  there  is  a  progressive  dis¬ 
tortion  of  these  areas  on  the  net  from  the  centre  outwards.  Fully  accurate  con¬ 
touring  requires  continual  changes  in  the  shape  of  the  counter  (Flinn,  1958).  One 
way  this  tedious  procedure  could  be  avoided  has  been  suggested  independently  by 
Noble  and  Eberly  (1964)  and  Muecke  (1965).  They  used  computer  techniques  to 
count  out  the  points  on  the  original  sphere.  A  grid  of  333  points  is  specified  in 
Muecke's  program  (by  their  direction  cosines)  and  when  a  reading  (also  given  by 
its  direction  cosines)  falls  within  a  cone  with  n%  cross  sectional  area  (n  =  1  -  9 
at  the  programmer's  choice)  around  a  grid  point,  the  counter  at  that  grid  point 
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is  given  an  increment.  The  values  of  the  counters  are  then  plotted  out  on 
an  equal-area  net. 

Because  the  computer  technique  elimates  plotting  errors,  only  errors  in 
field  observations  remain.  It  has  previously  been  shown  that  the  distribution  of  the 
probability  of  finding  the  true  pole  to  a  bedding  plane  might  reasonably  be  assumed 
to  have  a  Fisher  distribution  about  the  best  estimate  of  the  pole.  In  other  words, 
the  probability  of  the  true  pole  lying  at  angle  greater  than  0  with  the  best  estimate 
of  the  pole,  falls  off  with  increasing  0  according  to  the  Fisher  probability  density 
function,  P  (0) . 

The  probability  of  a  pole  with  a  Fisher  distribution  is  continuous  over  the 
hemisphere  of  projection.  Further  poles  plotted  will  increase  the  probability  of  a 
pole  at  any  locality.  This  will  be  the  sum  of  the  probabilities  of  individual  poles 
at  the  locality  (assuming  the  poles  are  disconnected  events).  For  instance  the 
probability  of  finding  at  least  one  pole  within  M  in  Fig.  4  is  the  sum  of  the 
probabilities  of  finding  A  there  and  B  there. 

Treating  probability  as  a  continuous  function  should  lead  to  some  reconsid- 
eration  of  techniques  of  contouring  it.  For  as  it  makes  little  sense  to  contour  a 
topographic  map  by  using  stations  determined  by  averaging  elevations  over  an  area, 
so  it  is  illogical  to  contour  probability  by  averaging  it  over  areas  as  large  as  1%  of 
the  hemisphere. 

Consider  the  probability  P  of  a  pole  in  a  small  area  JA  lying  between 


0  and  0O  from  P. 

2  /S 


/ 


where  P  is  the  probability  density  function. 

The  integration  is  over  the  surface  element,  6, A 


If  A  is  very  small  the  probability  has  a  constant  value  over  that  area. 
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Watson  and  Irving  (1 957)  show  that  i  f  K,  the  precision  parameter,  is  large, 
the  probability,  ,  of  the  true  mean  of  a  reading  lying  between  9^  ,  9^  of  the  esti¬ 


mated  mean  is  given  by 

{  -K(l  - 


X 


-e 


t 


-k(l  ~  xo 


The  probability  P  that  it  lies  in  the  area  £/\  is  given  by 


M  .  x' 

c‘/r  ( F-l l) 


if  9^ ,  92  are  small  and  the  hemis¬ 
phere  has  unit  radius,  and  c  =  1° 


in  radians 


Let  us  now  define  the  density  of  the  probability,  D,  as  the  probability 


per  unit  area  i 


X 


c2ty  (££ 

where  D  is  the  density  of  probability  in  the  small  area  between  9^ ,  and  92«  A  table 
(Table  5:1)  of  values  of  D  for  different  values  of  K  has  been  computed.  This  table 
also  includes  a  rough  weighting  of  the  probability  of  finding  the  reading  in  a  small 
area  at  differing  angles  from  its  estimated  position.  The  table  allows  the  adaptation 
of  Yanda  and  Muecke's  program  (Muecke,  1965)  to  increment  the  grid  of  counting 
locations  by  the  approximate  density  of  probaility  of  the  reading  in  the  small  area 

around  the  counter.  The  resulting  diagram  can  be  directly  contoured  to  appreciate 
the  distribution  of  probability. 

Whereas  the  relative  densities  of  probability  are  sufficient  to  show  the  dis¬ 
tribution  of  probability  on  the  diagram,  it  is  often  required  to  predict  the  number 
of  points  from  a  larger  sample  of  the  same  population  which  would  fall  within  a 
given  area.  This  is  the  object  of  constructing  contours  of  percentages  of  points 
per  percentage  area.  In  fact,  we  seek  a  number  N  such  that  N  x  R  (where  R  is 
the  percentage  of  the  diagram  within  the  N  contour  line)  is  the  number  of  points 
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enclosed  within  the  contour  line.  The  total  density  of  probability,  S,  associated 
with  1  point  is  given  by 


(5 


S  -  D  d  & 

> 0 

This  integral  can  be  approximated  by 

i=6 


where  9  is  the  angular  distance  at 
which  the  densi  ty  becomes  negligible. 


5 


D. 


I  -  0 


where  i  takes  integer  values 


is  calculated  in  the  table  for  various  values  of  K. 

Suppose  the  counters  are  so  placed  that  only  q%  of  the  density  of  probab¬ 
ility  is  recorded.  Then  the  value  of  the  density  at  each  counter  must  be  increased 

100  _ 2 

by  -  .  Throughout  the  table  a  factor  \\C  has  been  dropped  for  computational 

M 

convenience  so  that  if  D  is  the  value  of  a  particular  counter 

x 


2 


-  the  number  of  readings  falling 


within  the  R%  of  the  diagram  which  the  D  contour  encloses.  Yanda  and  Mueck's 

x 

program  has  therefore  been  redesigned  to  plot  out  a  second  diagram  giving  values 


of 


2  for  each  counter. 


The  procedure  outlined  has  some  analogy  with  that  suggested  by  Me  1 1  is 
and  advocated  by  Flinn  and  Friedmann.  They  suggest  drawing  a  circle  1%  of  the 
area  of  the  hemisphere  around  each  plotted  pole  and  contouring  areas  of  1 ,  2,  n 
intersecting  circles.  This  procedure  has  a  sound  basis  but  is  unable  to  distinguish 
between  the  extreme  cases  illustrated  below,  (Fig.  4).  The  extra  sensitivity  of 
the  computer  method  lies  in  its  ability  to  distinguish  between  cases  like  this,  its 
elimination  of  plotting  errors,  and  its  sensitivity  to  varying  K  values  of  the  input. 
Fig.  5  shows  this  schematically. 
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TABLE  5:1  VARIATION  OF  DENSITY  OF  PROBABILITY  WITH  ANGULAR 
DISTANCE  FROM  THE  ESTIMATED  MEAN  OF  THE  SAMPLE 
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